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We present Bg ^ p, Bg ^ co, Bq ^ K*, Bg —t K* and Bg ^ (j) form factors 
from light-cone sum rules (LCSR) at 0(ag} for twist-2 and 3 and 0{a^) for twist-4 
with updated hadronic input parameters. Three asymptotic light-cone distribution 
amplitudes of twist-4 (and 5) are determined, necessary for the form factors to 
obey the equations of motion. It is argued that the latter constrain the uncertainty 
of tensor-to-vector form factor ratios thereby improving the prediction of zeros of 
helicity amplitudes of major importance for B —)> K*£i angular observables. We 
provide easy-to-use hts to the LCSR results, including the full error correlation 
matrix, in all modes at low as well as combined hts to LCSR and lattice results 
covering the entire kinematic range for Bg —?• 77*, Bg —?• K* and Bg —)• (/>. The 
error correlation matrix avoids the problem of overestimating the uncertainty in 
phenomenological applications. Using the new form factors and recent computa¬ 
tions of non-factorisable contributions we provide Standard Model predictions for 
B —)• 77*7 8-8 well as 77 —?• 77*7"’“7“ and Bg —)• 4>p~^pL~ at low dilepton invariant 
mass. Employing our B —?■ {p,u}) form factor results we extract the CKM element 
|Uub| from the semileptonic decays B —)> {p,uj)£iy and hnd good agreement with 
other exclusive determinations. 
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1. Introduction 


Exclusive semi-leptonic B decays are important tools to test the Standard Model (SM) and 
to look for new physics. Among these processes, the decays B —)• K*{-^ it'7r)/x+/r“ and 
Bg —>• (/>(—)• K~^ij,~ are of particular relevance as their angular distributions give access 
to a host of observables that are sensitive to new physics (e.g. [1] for arecent review). Predicting 
these observables, either within the SM or beyond, requires the knowledge of the form factors 
(FFs) “ in the case of S —)• F transitions, these are 7 functions of the dilepton invariant mass 
squared q^. In the low region, where the vector meson is energetic, the FFs can be computed 
using the method of sum rules on the light cone (LCSR) whereas at high q^ the FFs can be 
computed using lattice QCD. 

In this work we present an update of the FF computation in [2], for the modes Bq —)• p, 
Bq Lo, Bq ^ K*, Bg — )■ K* and Bg ^ (p (with q = u,d), using current hadronic input and a 
concise discussion of the role of the equation of motion (EOM) in correlating vector and tensor 
FFs. The FFs are fitted to the z-expansion parameterisation in the helicity basis, retaining all 
correlations among the expansion coefficients^. This information is made publicly available as 
ancillary files on the arXiv web pages in a form which is easy to use for phenomenology. 

Crucially the correlation of the uncertainties avoids overestimating uncertainties in observ¬ 
ables. A particularly important example are the angular observables in i? —)• K*p~-type, 
decays since they are sensitive to ratios of FFs and zeros of helicity amplitudes. For the latter 
two, the uncertainty is considerably reduced when taking correlation into account. 

We argue, extending the work in Ref. [4], that the use of the EOM enforces the correlation 
of the non-parametric, sum rule specific, input parameters. This can be seen as an application 
of the large energy limit (LET) ideas [5] to the sum rules on the light-cone. It is in giving 
numerical predictions and not relying on the heavy quark limit that the LCSR computations 
go beyond the LEL ideas [5] (this includes the case factorisable hard Q;<j-corrections [6]). The 
LCSR therefore give corrections to the LEL [5] and soft-collinear effective theory (SCET) [6] 
relations. Going beyond the SCET framework of two soft FFs and hard a^-correction in 
the heavy quark limit involves using the numerical predictions from LCSR, e.g. [7,8]. Going 
beyond the SGET framework has become increasingly important since observables designed to 
minimise the impact of the soft FFs [9] are, of course, sensitive to l/rub-corrections. 

We perform combined fits of the FFs to the LGSR at low and a recent lattice computation 
at high (fi [10,11]. This serves to test the consistency of the two complementary methods and 
provides FF sets valid over the entire kinematical region. We extract the CKM element JIFb] 
from B —)■ (w, p)iv BaBar- and Belle-data using the B —)■ (w, p) FF predictions of this paper. 
This can either be viewed as an extraction of iFjb] or as a check of the normalisation of the 
FF when compared to global fits or R —>• Triu extractions of iKibl- In addition to the FFs, 
the calculation of R —)> Vi^i~ observables involves non-factorisable contributions from the 
weak hadronic Hamiltonian. Some of these contributions have been recently computed within 
LCSR. Including all of these ingredients, we present SM predictions for the branching ratios 
and angnlar observables of R —)■ K*p~ and Bg —)• 4)p'^p~. We also compare the prediction 
for the branching ratio of R —iL* 7 , that has been measured precisely at the R factories, to 
the data. 

The paper is organised as follows. In section 2 we present and discuss the seven R —?■ R FFs 
within the LCSR context, discussing the implication of the EOM , finite width effects, input 


^Similar fits retaining correlations have recently been performed for the B —>■ tt FFs [3]. 
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parameters and the interpolating fits to the lattice data. In section 3 phenomenological aspects 
of i? —)• B —)■ K*'y, Bg —)• B —)• K *versus Bg —)• and the determination of 

iVubl from B —)• {p,uj)(.v are discussed: see subsections 3.1, 3.2, 3.3, 3.4 and 3.5 respectively. 
Conclusions figure in section 4. Appendix A assembles aspects of the EOM, explicit tree level 
results, scheme dependence and remarks on the Borel parameters. The determination of the 
three light-cone distribution amplitudes Ay (twist-4) and (twist-5), using an alternative 
method, is discussed in appendix B. A detailed discussion on the determination of the decay 
constants from experiment is given in appendix C. Conversion between bases, further plots 
and fit coefficients are given in appendixes D, E, and E respectively. The effect of the sizeable 
B^-lifetime is worked out in appendix G. 


2. B ^ V form factors from light-cone sum rules 


The short distance matrix elements, relevant for the dimension six effective Hamiltonian, are 
parameterised by seven EFs^ 


(K*(p,7/)|s7^(1T75)6|S(pb)) = P>^Vi{q^)±P^V2{q^)±P^V^{q^)±P^Vp{q^) , 

{K*{p,p)\-siq,a^^^{l ± l^)b\B{pB)) = P'^T^{q^) ± P^T 2 {q^) ± P^T^{q^) , (1) 

where the Lorentz structures Pj^ are defined as in [12] 

= i{v* • q)q>^ , Pf , (2) 

P 2 = ■Q)ip + PbT} , P 3 =i{v*-q){q'" -^— ^{p + PbV} , 


with the eoi 23 = +1 convention for the Levi-Civita tensor. The relation Pi(0) = P 2 ( 0 ) holds 
algebraically. The parameterisation (1) makes the correspondence between vector and tensor 
EEs explicit. The correspondence of the to the more traditional EEs ^ 07 , 2,3 and V is 

as follows 


Vp(q^) = ^^^* Ao(q^) , Vi(q^) = 


M ! 2\ I'^B + rriK* , ( 2\ 

V3(g ) = (■- 2 -) - 


-V{q^) 
mp + rriK* 
mp — rriK* 


V 2 {q^) = 


mp — mx* 


q- 




q 

p — 


(3) 


The relation ^ 3 ( 0 ) = Ao(0) assures finite matrix elements at q^ = 0. The last relation in (3) 
indicates that one FF out of A 177 is redundant.^ The pseudoscalar matrix element is related 
to Aq through an axial Ward Identity: 




(4) 


The projection on the helicity basis, using the Jacob-Wick polarisation convention, is given 
in appendix D. In the next section we briefly discuss the use of the method of LCSR before 
investigating the implications of the EOM on certain sum rule specific parameters. 

^Due to the composition of the wave functions of the p® ~ l/y/2{uu — dd) and lu ~ l/\/2{im — dd), extra 
factors have to attached to the matrix elements on the left-hand side, cf. [2]. They read: c“o = = 

cZ = ct = V 2 and cv = 1 in all other cases. 

®From the viewpoint of the projections the traditional nomenclature is unfortunate. It would have been better 
not to have A 2 at all and use the notation Ai —>■ A 2 . 
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2.1. Calculation of the form factors in light-cone sum rules 

Light-cone sum rules (similar to QCD sum rules [13,14]) for the FFs are derived by considering 
the correlator of the time-ordered product of two quark currents evaluated between the final 
state on-shell meson (in this case V) and the vacuum [15,16]. On expanding this correlator 
about the light-cone, one obtains a series of perturbatively calculable hard scattering kernels 
convoluted with non-perturbative, universal light-cone distribution amplitudes, ordered by in¬ 
creasing twist (dimension minus spin). Reasonable convergence of the LC-expansion is formally 
and by experience limited up to ~ ©(mbAgcu) — 14 GeV^. In the hadronic picture the cor¬ 
relator is expressed as the sum over excited states, the dominant state being the R-meson, and 
this is followed by the continuum. Assuming quark-hadron duality above a certain continuum 
threshold sq [13,17], an approximation referred to as the semi-global quark-hadron duality 
assumption, one arrives at an expression for the lowest lying hadronic parameter in terms of 
an expression of partonic QCD and sq- A Borel transformation which ameliorates both the 
hadron and the parton evaluation leads to a numerical improvement of the procedure. 

Light-cone sum rules results, with light-meson distribution amplitudes (DAs), have been 
computed for the B ^ P transition up to twist-3 0(as) in [18-20] and for the B ^ V 
transition up to twist-4 at tree level and twist-2 0(as) [21] as well as twist-3 0{as) [2]. In 
this paper we make use of the results in [2].^ Alternatively the FFs can be determined using 
R-meson DA and an interpolating current for the R-meson. This program has been pursued 
in [22] at tree level in QCD and in SCET [23]. In this work we improve on the previous LCSR 
work [2] by 

• computation of the full twist-4 (and partial twist-5) 2-particle DAs contribution to the 
FF (appendix B - available a downloadable Mathematica notebook), 

• determination of the DAs (twist-5), in the asymptotic limit, filling a gap in the 
literature (appendix B), 

• discussing the impact of the EOM on uncertainty correlations (section 2.2), inclnding the 
aspect of scheme-dependence (appendix A.5), 

• explicit verification of the EOM at tree level (appendices A.3,A.4) , for the asymptotic 
2-particle DAs including 0(ms)-corrections. 

• verification of the compatibility of the composite operator renormalisation with the EOM 
(appendix A.5.1), 

• discussion of non-resonant background for vector meson final states (section 2.3), 

• determination and usage of updated hadronic parameters (section 2.4), specifically the 
decay constants (appendix C), 

• fits with full error correlation matrix for the z-expansion coefficients (section 2.5), as well 
as an interpolation to the most recent lattice computation (section 2.6). 


^In [2] the size of the twist-3 0{as) corrections were not explicitly given. At = 0 the twist-3 0{as) 
corrections lead to a raise of around 10% of the FF Ti. 
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2.2. Equation of motion and form factors 

In this section we reiterate the use of the EOM [4]. As discussed in [4] this is of importance in 
reducing the uncertainty between certain FFs. Below we give more details and strengthen the 


argument. The following FOM 

^— 

id'"{siaf,^{j5)b)= - {ms±mb)sj^{'y5)b + idij,{s{j5)b)-2siD^j, ( 75 )^, (5) 

are valid on physical states. Equations of the form (5) are sometimes also referred to as Ward 
identities. In particular, evaluated on (F| .. .\B), Eq. (5) yields 

Ti{q‘^) + {nib + ms)Vi{q^) +T>i(g^) = 0 , (6) 

T2{q^) + {nib - ms)V2{q‘^) + ^> 2 ( 9 ^) = 0 , (7) 

Tsiq^) + {nib - ms)V3{q‘^) + V^{q^) = 0 , (8) 

{mb - ms)Vp{q‘^) + (vp{q^) - - -Vp(g^)^ = 0 . (9) 

V mb + ms ) 


One of the above four equations corresponds to each of the directions ( 2 ) [4,24], where Pj are 
defined 

{K*{p,T])\-s{ 2 ibr{\±-ib)h\B{vB)) = P^Vr{q^)±Pl^V 2 {q^)±P!^V^{q^)±Pp 2 p{q^) , (10) 

in complete analogy with (1). Note that the id^(s( 75 ) 6 ) operator only contributes to Pp ~ 
since the total derivative is replaced by the momentum transfer q^. In Eq. (9) we have included 
this contribution into round brackets with the other derivative FF. Before discussing Eqs. ( 6 )- 
(9) in various limits, we wish to stress that the equations are completely general and have to 
be obeyed by any FF determination. 

The Isgur-Wise relations [25] follow from a clear physical picture. At low recoil the non- 
relativistic heavy quark effective theory applies and it can be shown that Pj are suppressed by 
(AgcD/^TTfe) with respect to the vector and tensor FFs [24]. This raises the question of whether 
there are combinations of Pj’s which are small at large recoil. Eqs. (6,7) are direct candidates 
but Eqs. (8,9) require some thought because of the common direction q^,. In fact in Eqs. (8,9) 
the poles in q^ cancel between the FF VsjVp and P^,Pp which implies that V^^Vp are not 
individually small. Since the (AT* [ 57 ^ 7561 . 6 ) matrix element is free from singularities, adding 
Eqs. (8,9) yields a combination for which the derivative FFs are potentially small. We define 

_ {mb + ms)Vi{q^) _ mb + ms V{q^) 

Ti{q^) mp + mK* Ti{q^) ' 

_ {mb - ms)V 2 {q‘^) _ mb - ms Ai{q^) 

T 2 {q^) mB-mK*T 2 {q^) 

_ {mb - ms){V 2 {q^) - C23(g^)(V3(g^) + Vp(g^))) 

T2{q^) - C23{q^)T3{q^) 

- - ms){Vo{q‘^) - C23{q^)Vp{q^)) 

ro(V) 

np - mx* T 2 {q^) - C23{q^)T3{q^) 


the following ratios 
r^{q^) = 

r\\{q^) = 

ro+t{q^) = 


6 












Figure 1: Plots of rj_, ry and ro+t Eq. (11) (for the B —>• i^*-transition) as a function of . The 
deviation from unity measures the relative size of the derivative form factor with respect 
to the tensor and vector FFs. The mass used in (11) is the pole mass, an issue discussed 
in appendix A.5. For the explicit mb mass we use the central value and do not include 
an error since in the B —>■ V£i helicity amplitudes the Wb-pole is not present. The fact 
that the quantities tend towards one for very high is expected from the viewpoint of the 
Isgur-Wise relations and proves a certain robustness of the LCSR-results for high q^. Similar 
ratios have been plotted in previous work in the context of Isgur-Wise relations [21] and LET 
relations [68] but the observation that this might be useful for the correlating the continuum 
thresholds was not made. In particular the derivative form factor were not identified as a 
independent objects. 


where Xq = X 2 — C 23 X 3 for X = T,V with 023 ( 5 ^) being a kinematic function defined in(D.3). 
The deviations of these quantities from one measure the relative size of the derivative FFs with 
respect to vector and tensor FFs, 

_ 1 _L — 1 _L ^11 — 1 _L not 

r_L — 1 + , ry — 1 -b — , ro+t — 1 H — , ( 12 ) 

J-l £2 2 o 

where 'D± = Di, Py = 1)2 and T>o+t = 1^2 — C23{'1^3 + T’p). In Fig. 1 we show plots of these 
ratios from 0 < < 14GeV^. The quantities r+y and, somewhat less, ro+t are found to be 

very close to unity over this range. The basic idea is that if the 'Di are considered as regular 
FFs with controlled uncertainty^ then this implies a high degree of correlation between vector 
and tensor FFs of a given polarisation. This is partly reflected in the controlled error bands. 

The aspect of the correlation between the continuum thresholds is discussed in some more 
detail in appendix A.2. Here we just summarise the main argument and result. Based on the 
EOM we argue, conservatively, that the continuum thresholds for tensor and vector FF cannot 

®For the B —> K* channel at = 1 GeV^ the corrections due to twist-4 and Ua-correction for 
{{Ti,T>±), (To,Vo+t)} are {(4, 6), (7, 28)}% and {(12, 27), (11, 31)}% respectively indicating regularity of Vi 
with regard to the twist- and the a^- expansion. 
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differ by more than 1 GeV^ in order for the continuum thresholds of the derivative FF not to 
take on absurdly low or high values. This argument is less compelling for the (0 + t)-direction, 
as can be inferred from the plots, resulting in lower correlation and larger error bands. We 
stress that if we were to impose standard error bands, say sq = 35(2) GeV^ for the sum rule of 
the T>i FFs then the error bands for rj-ratios would shrink to the 1%-level. 

The smallness of the derivative form factors (cf. also [4] for further references and more 
physical discussion) is related to the findings of Gharles et al [5] within the LET and its 
extension into SGET [6,26]. The similarity is the use of the EOM of motion but the difference 
is that in this work the EOM are directly implemented within QGD whereas in the prior 
work the EOM are used at the level of an effective theory in l/m;,. This results in differences 
at the level of l/rrib corrections. For example in [6] the ratio analogous to r±_ (11), which 
we shall denote by = mB/{mB + rriK*) iy/Ti) = 1 + C>(m^^, asm°), differs from ours 
by niB —>■ {mb + mg) which is indeed 0{m'j^^). For completeness let us mention that the 
symmetry breaking corrections to were computed to 0{asm^) and 0{a‘lm^) in [6] and 
[27,28] respectively. Even though the twist expansion of LGSR is not a l/rub-expansion, as 
first stated for B ^ P and B ^ V decays in [16] and [29], it is of interest to examine the various 
twist-quantities from the viewpoint of the standard heavy quark scaling prescriptions [16]. At 
the level of all explicit calculations in the literature it is found that ri( 0 ) ~ E(0)|twist-2,3 ~ 

and ri(0) ~ E(0)| twist-4 ~ m^ The derivative EFs follow Pi (0)| twist-2,3 ~ m^ [4] 
which is in agreement with the explicit computation in [5]. For this work we have explicitly 

_■? /o _Q /o 

checked that Pi( 0 )|twist -4 ~ m^^ ' and that Pi(0) ~ 0{as)m^^ ' , in accord with the results 
of Beneke and Feldmann [6]. In summary we may state that the parametric statements and 
the previous numeric statements give a consistent picture. 

2.3. Discussion of non-resonant background effects 

The signal final state in a P —)• p{-^ 'K'K)lv-iype decay, serving as a template for any B —)• Vli (.2 
decay, is tt'kIv. Hence in principle the decay ought to be analysed via a H —)• vrTr type form 
factor.® The analysis of H —)• ptv therefore becomes a matter of how background, finite 
width and S,P-wave effects are treated or discerned. This question arises in any theoretical 
computation as well as in any experimental measurement. It is therefore important that both 
theory and experiment treat these issues in a consistent way. 

Let us contrast the p hnal state with the Trvr-state from a pragmatic viewpoint relevant for 
this paper. The orbital angular momentum of the vr-Tr-state is either / = 0,1,.. (S', P, ..-wave). 
If the vrvr-state originates from a p-meson then it is necessarily in a P-wave state and shows 
a distinct angular distribution. Hence this contribution can be separated through an angular 
analysis from other type of partial waves. ^ We therefore conclude that the S-wave contribution 
is not to be included in a P —)• p FF and we therefore do not need to attribute any additional 
uncertainty to it. 

We turn to the question of the treatment of the P-wave. For the sake of concreteness we 
discuss the B —)• J/Ttf+tt" measurement of the LHG collaboration [37]. In a certain window 
around ~ mj, the Trvr-spectrum, in the P-wave, is fitted through an ansatz of a resonant p 

®Within the framework of LCSR this could be done by using the a two-pion DA [30-33]. The technology for 
pursuing a i? —>■ tttt FF computation on the lattice has been put forward recently in reference [34]. 

^The importance of separating the 5-wave, in the context of B —^ K*££-type decays, was emphasised not 
long ago in [35]. Thereafter the 5-wave FFs B —> (Avrls-wave were computed in LCSR in the tree level 
approximation [36]. 
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and the two excited p(1450) and /?(1700) states. The main fit parameters are the amplitudes 
(complex residues) of the Breit-Wigner ansatz whose values determine the interference pattern. 
A non-resonant P-wave background is usually not fitted for since it is assumed that the S-wave 
is dominant in the non-resonant background.® This raises the question of how a theoretical 
framework like LCSR can accommodate this complex procedure. The answer is surprisingly 
pragramatic. As long as the experimental input into the LCSR is treated in the same way 
there is no systematic effect. 

Let us argue this point in some more detail. In current LCSR determinations of the B ^ V 
FFs the F-meson is described by vector meson DAs. The latter are mainly characterised by 
the longitudinal and transverse decay constants fy and fy respectively and it is therefore 
important to know how they are obtained. The method of choice for determining fy is exper¬ 
iment: e'''e“ —)• F°(—)• PP) (for = p^,(j),uj) and —)• F'''(—)• PP)v (for F+ = 77*’+, p"*") 
cf. appendix C. As long as the experimental treatment of the p-meson versus vrvr-signal event 
is the same as for the semileptonic decay P —)• (p —?■ Tnr)iL' the decay constant encodes the 
same definition of the p-meson as used in P —)• pii'. The other quantities associated with the 
p-DA are not directly accessible in experiment. For example the transverse decay constant fy 
is obtained through ratios of fy/fy from lattice QCD and sum rules where one would expect 
effects of the treatment of the p-meson to cancel to a large extent or to be taken care of by the 
respective uncertainties. 

We may discuss the same reasoning from the viewpoint of a computation using a two-pion 
DA instead of a p DA. Let us consider for example the contribution of the DA that couples 
to the vector current. The latter is described at lowest order in the conformal spin expansion 
by /]j times the asymptotic DA which follows from conformal symmetry. From the formulae 
in [31] it is seen that the analogue of the two-pion DA is given by the pion FF f^\q‘^) times 
the asymptotic DA. Somewhat symbolically the transitions in terms of p and 2-pion DAs are 
given by 


vector p—DA : 

/II • BW(mL) 

^ ^ fp'lTT’pdpmr 1 

ml^-mj-impTp 

(13) 

vector TTTT—DA : 

f+\ml^) 

^ fpl^pdpTTTV 

^ 9 0 • T-' ’ 

-m^p- irupYp 

(14) 


for ~ m?p where gp-j^w is the p —)• tt-tt decay constant, BW stands for some type of Breit- 
Wigner ansatz and the dots stand for all contributions other than the p-resonance from the tt-tt 
P-wave. Our argument is that unless one is specifically interested in the local m^^-behaviour 
this contribution can and is effectively absorbed into fp upon integration over the p mass 
window in the experimental analysis. For higher order conformal spin corrections, i.e. higher 
Gegenbauer moments, and other decay constants the same reasoning applies. The strong 
rescattering phases in the TTTT-channel are universal in each partial wave and do not distort the 
result. 

In summary, from a pragmatic viewpoint as long as the treatment of the p(—)• 7r7T)-meson 
is the same that is used for the extraction of fp, the LCSR should not suffer from sizeable 

®In the cases where a background has been searched for in i? —> p£u, it has been found to be consistent with 
zero [38,39]. Whether or not future experiments can discern the background is difficult for us to judge but 
we argue that from a pragmatic point of view this might not be necessary. 
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/"[GeV] 

/^[GeV] 

a| 


a| af 

cl [46] 

p 

0.213(5) 

0.160(7) 

0.17(7) 

0.14(6) 

- - 

0.030(10) 

Ijj 

0.197(8) 

0.148(13) 

0.15(12) 

0.14(12) 

- - 

idem 

K* 

0.204(7) 

0.159(6) 

0.16(9) 

0 .10(8) 

0.06(4) 0.04(3) 

0.023(8) 


0.233(4) 

0.191(4) 

0.23(8) 

0.14(7) 

- - 

0.024(8) 


Table 1: The determination of is discussed in some detail in appendix C. The fine structure constant 
a, relevant for the extraction of fy, is evaluated at /r = 1 GeV ~ my- Scale dependent 
quantities, e.g. and C|j are evaluated at = 1 GeV. The parameters are 

taken to be the same as in [47] which include computations from [46,48-51]. The decay 
constants are obtained from /H through ratios r[X] = /j^(2 GeV)//]|- with r[p] = 0.687(27), 
r[K*] = 0.712(12) and r[(j)\ = 0.750(8) from lattice QCD [52]. For the w-meson we adopt 
r[a;] ~ r[p] in view of a lack of a lattice QGD determination of this quantity. Twist-3 DA 
parameters are taken from the values for , Wg and cDg [46] which include cl (quoted in the 
table), a;|, . The twist-4 3-particle DA parameters are neglected since they are at the sub 

per mill domain. Again, for the w-meson we adapt the same values as for the p-meson since 
a specific determination is lacking. 


additional uncertainties. It therefore seems that in practice the uncertainty is a small frac¬ 
tion of the the P-wave background which itself is around 5%.® In view of all other sizeable 
uncertainties we refrain from adding any further error due to this effect and reemphasise the 
importance of comparing our result only with the P-wave contribution of the corresponding 
TTvr-pair. Whereas the analysis in this section questions the practical impact of using a two-pion 
DA around the the p-meson mass, it is of course interesting to look at the decay B —?• Tmlv 

in other regions of phase space. For recent theoretical developments of B^^ we refer the reader 
to [43,44] which are though not yet at the level of maturity of [45]. 

2.4. Input parameters and uncertainties 

The uncertainty of the LCSR results for the FFs is determined from the uncertainties on the 
input parameters, the factorisation scale p,F and the Borel parameter as well as the effective 
continuum threshold sq- The values of input parameters used in our calculation, along with 
the errors assigned can be found in Tab. 1. We draw the reader’s attention to the fact that it is 
the quantity [F{q‘^) ■ fs], where F stands for any of the seven form factors, that is determined 
from the correlation function. Therefore one needs to divide by fs in order to obtain the FF 
P. It is well known and appreciated that the uncertainty in as is considerably reduced when 


^Despite this aspect it is of interest to estimate the non-p background. One can get an idea by analysing 
the pion FF f^\q‘^){pi — P 2 )/j = (^(Pr)7r(p2)|i[l™|0). A measure of the non-resonant background around 
the p-meson peak is given by the difference of the model-independent determination of the pion FF using 
data on TTTr-scattering phase shifts and the Omnes-dispersion relation versus a fitted p-meson Breit-Wigner 
ansatz. Around the p-meson mass window the difference is found to be 5% [40]. Similar conclusions can be 
reached when considering the figures in [41] with and without the p(1450) and p(1700) contributions. We 
note that 5% is of the same order as the S-wave background found in i? —>■ Tmiv [42]. 
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sum rule in is taken to same order as for the quantity [F ■ f^] ■ For example fs increases by 
~ 9% at 0{a^gj5Q) whereas the combination (/+^^/b)lcsr/(/b)sr only increases by 2% [53]. 
Therefore we make use of the QCD sum rules result at 0(as) [54,55] for fs- 

The two sum rule specific parameters are the Borel parameter and the effective con¬ 
tinuum threshold sq- For reasons of consistency the Borel parameter is to be chosen at an 
extremum (cf. appendix A.6) which serves as a quality control parameter. The continuum 
threshold is more problematic and the final result does depend on the choice. Hence our recipe 
for the error analysis is to assume a sizeable uncertainty for the continuum threshold. The 
new ingredient of our analysis is that we have argued that the EOM results in correlations 
between continuum thresholds of certain FFs; (cf. appendix A.2 for an elaborate discussion). 
The correlations used are summarised in and in between Eqs. (A.8) and (A.9). The correlation 
of the continuum thresholds are such that the relative uncertainty is 1 GeV^ which has to be 
compared to the individual uncertainty of 2 GeV^ or the uncertainty of their sum which is close 
to 4 GeV^. The influence of the Borel parameter on the light-cone sum rule is negligible as 
compared to the continuum threshold of the light-cone sum rule and we therefore do not vary 
them separately for each FE. The Borel parameter dependence of the fs-sum rule is sizeable 
and is taken into account and contributes to the uncertainty of the normalisation of the FEs. 
The intermediate states for the light-cone and the /b sum rule are the same since they couple 
to the same interpolating current Jb- It would therefore seem absurd, or contradictive to 
the method, if the corresponding continuum thresholds were far apart. We implement this 
reasoning by correlating Sq® and at the 50%-level which implies that the uncertainty on 
the difference is 2GeV^; a factor of \/2 lower than without correlations. 

We turn to the choice of the actual central values of the continuum threshold and the Borel 
parameter. It is useful to recall that if the sum rules were perfect then the LCSR FF would be 
independent of the Borel parameter. In reality a small Borel parameter is desirable from the 
viewpoint of suppressing any higher states in the spectrum whereas a large Borel parameter 
improves the convergence of the light-cone operator product expansion (LC-OPE). In practice 
one therefore chooses a compromise value which is usually found as an extremum. The flatness 
of the FF around this extremum as a function of the Borel parameter is a measure of the quality 
of the sum rule. In appendix A.6 it is shown that extremising in the Borel parameter is formally 
equivalent to imposing a daughter sum rule for m^. From the viewpoint of the physics, the 
effective continuum threshold is expected to he somewhere between (ms -|- 2m7r) — 30.9 GeV^ 
and {mB + 'irip)'^ ~ 36.6 GeV^ with the true value being closer to the latter since the production 
of a p-meson is much more likely than the production of two non-resonant pions. The twist-4 
contribution for ±-helicity (Tip, V, and Ai) is around 5% whereas for 0-helicity {T 21 ,, A 12 and 
Aq) they are just below the 10%-range. Guided by the relative size of the twist-2 and twist- 
3 radiative versus tree contribution^*^ we estimate the uncertainty due to the missing 0(as) 
twist-4 contribution^^ by associating a Gaussian error of 50% to the latter. 

In order to limit contamination due to higher states we verify that the continuum contribu¬ 
tion does not exceed 30%. If one assumes that semi-global quark hadron-duality itself works 
at the 30%-level the additional suppression reduces this error to just below the 10%-level. The 
sum rule parameters, with some more details in the caption, are given in Tab. 2. 


remind the reader that the actual impact of the radiative corrections on the FF result is considerably 
smaller since a large part is absorbed by the radiative corrections to fs (cf. the beginning of this subsection). 
^^More precisely no 0{my) are included at 0{as). (cf. table II in [2]). We impose a 50% uncertainty on the 
corresponding tree-level terms. 
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Bq 

JBq 

]\j2 

„LC 

Bd 

4.1(1) 

34.2(2) Cc/{u)g2Mj^^ 

35(2) 

Bs 

4.4(1) 

35.4(2) Cc/{u)g2Mj^^ 

36(2) 


Table 2: Sum rule parameters for and Bg sum rules. All numbers are in units of GeV^, 

and Sq ' denote the Borel parameter and continuum threshold of the fs sum rule and 
the LCSR of (where F stands for a FF) respectively. The difference between the 

Bd and Bg continuum thresholds follows {msg + A)^ = sols^ and + A)^ = SqIb,■ The 

average momentum fraction of the transition quark {u)q 2 (cf. [19] for the definition) varies 
smoothly from 0.86 at = OGeV^ to 0.77 at = 14GeV^. Dividing the sum rule parameter 
by this quantity serves to take into account g^-dependence the Borel parameter under the 
extremisation procedure. The value Cc = 2.2 is determined through the mentioned procedure 
of extremisation. The criteria in the text imply that the Borel parameter of the LCSR is 
considerably higher than that from the /^-sum rule [19]. 


Fi 

JP 


m^ 7 /GeV 

^0 

0 " 

5.279 

5.366 

Ti,R 

1 “ 

5.325 

5.415 

F 2 , T 23 , Ai,Ai2 

1 + 

5.724 

5.829 


Table 3: 


Masses of resonances of quantum numbers as indicated necessary for the param- 
eterisation of FF F) for b ^ d and 6 —)■ s transitions. 


2.5. Series expansion fits to LCSR form factors 


As mentioned in the introduction, for phenomenological analyses of rare decays, it is crucial to 
take into account the theoretical uncertainties of the B ^ V FFs and the correlations among 
them. In order to facilitate the use of the LCSR results, we perform fits of the full analytical 
result to a simplified series expansion (SSE), which is based on a rapidly converging series in 
the parameter 

\/^+ “ ^ “b y/t+ ~ ^0 



where t± = (m^ ± mvY to = ^+(1 ~ 


We write the FFs as 


Fi{q^) = Pi{q^) “fc “ ^(0)] ^ ’ (16) 

k 


where Pi{(f‘) = (1 — is a simple pole corresponding to the first resonance in the 

spectrum. The appropriate resonance masses are given in Tab. 3. We consider fits that are 
truncated after the quadratic term in z, i.e. we will have three fit parameters 00 , 1,2 for each of 
the seven FFs. We will see in section 2.6 that a three-parameter fit is sufficient for a combined 
fit to lattice and LCSR results in the entire kinematic range relevant for B —)• decays. 

Note that the parameterisation (16) differs from that used in [10,11]. It has the advan¬ 
tage that the value of the FF at = 0 is among the fit parameters, F)(0) = Oq. We 
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B - 


K* 

B 


• P 

B 

—> 

UJ 

Bs 

- 

^ 

Bs 


K* 

Am 

0.356 

± 

0.046 

0.356 

± 

0.042 

0.328 

± 

0.048 

0.389 

± 

0.045 

0.314 

± 

0.048 

Am 

0.269 

± 

0.029 

0.262 

± 

0.026 

0.243 

± 

0.031 

0.296 

± 

0.027 

0.230 

± 

0.025 

^12(0) 

0.256 

± 

0.033 

0.297 

± 

0.035 

0.270 

± 

0.040 

0.246 

± 

0.029 

0.229 

± 

0.035 

F(0) 

0.341 

± 

0.036 

0.327 

± 

0.031 

0.304 

± 

0.038 

0.387 

± 

0.033 

0.296 

± 

0.030 

ri(0) 

0.282 

± 

0.031 

0.272 

± 

0.026 

0.251 

± 

0.031 

0.309 

± 

0.027 

0.239 

± 

0.024 

^2(0) 

0.282 

± 

0.031 

0.272 

± 

0.026 

0.251 

± 

0.031 

0.309 

± 

0.027 

0.239 

± 

0.024 

^ 23 ( 0 ) 

0.668 

± 

0.083 

0.747 

± 

0.076 

0.683 

± 

0.090 

0.676 

± 

0.071 

0.597 

± 

0.076 


Table 4: Values of the FFs at = 0 and their uncertainties. The tensor FFs are renormalised at the 
pole mass of the 6-quark /ruv = 4.8 GeV. For a more detailed error breakdown we refer the 
reader to the table 7 of the previous LCSR FF work [2] . 


prefer this parameterisation as it allows to impose the exact kinematical relations Ao( 0 ) = 
{SmBmv )/(which is equivalent to Ao( 0 ) = ^ 3 ( 0 )) and ri( 0 ) = T 2 ( 0 ) at the 
level of the SSE coefficients as 


a 


Ao _ 

0 ~ 


SmBmv 

2 T 

— niy 


a, 


Ai2 


a 




(17) 


The results for the FFs at = 0 are provided in Tab. 4. To determine the fit coefficients ctj, 
the uncertainties, and the correlations between them, we first generate an ensemble of = 
500 input parameter sets where the values of the input parameters are randomly distributed 
according to a multivariate normal distribution with the location given by the central values 
and the covariance given by the uncertainties and correlations of the input parameters discussed 
above. We then compute all FFs at integer values of between 0 and 14 GeV^. Finally, we 
fit the z expansion to all seven FFs for the N ensembles of FF values and extract the mean, 
variance, and correlation of the z expansion coefficients ao,i, 2 - Since we impose the exact 
conditions (17) throughout, the number of independent fit parameters is 19. 

The resulting mean and variance are shown in Tab. 14. We do not reproduce the full 21 x 21 
correlation matrices in the paper but rather provide them as downloadable ancillary files which 
are available on the arXiv preprint page (see appendix F for details). Here we merely note 
that these correlations are sizeable and it is crucial to include them when using the FFs in 
phenomenological analyses. 

With these results at hand, the uncertainty of an observable T (e.g. angular observable) can 
be computed via 


u2(T) = 


E 


d^{Fi) 

dal 


cov(al,aj) 


daj 


(18) 


where i,j = 1... 7 denotes the FF index and k,l = 0.. .2 parameterises the expansion coeffi¬ 
cients of the z-series. The covariance matrix is defined as 


cov{al, aj) = coTT(al,aj) a(al) a(aj) (no sums) 
in terms of the correlation matrix and the variances. 


(19) 
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As an example let us write the formula relevant to the ratio TZi{q^) = (raB^vny) jjV{q^) 
whose difference from 1 marks difference from the large energy limit [6]. At = 0 the error 
of the FF ratio is given by 


a 



( ^ 


2 



a 


Ti 

0 


f a{al) 

V 


- 2 - 


,17 a, 


--corr(aoi,ao ) 


“o “o 


( 20 ) 


from which we obtain = 0.97 ± 0.04. 


2.6. Interpolating between lattice and LCSR form factors 

The LCSR and lattice FF calculations are complementary since the former is valid at low (f' and 
the latter at high q^. Performing a combined fit of the SSE parameterisation to both lattice and 
LCSR results is useful for two reasons. First, whether a good fit to two completely independent 
methods in two different kinematical regions is possible at all is a powerful consistency check 
of those methods. Second, in phenomenological analyses constraining physics beyond the SM 
combining both observables at low and at high q^, one needs a consistent set of FFs for the 
full (f‘ range. 

To obtain this combined fit, we first generate pseudo-data points with correlated theoretical 
uncertainties at three q^ values both at low and at high q^. For LCSR at low q^, we proceed as 
in the previous subsection. For the lattice FFs at high we make use of the parameterisation 
of lattice FFs provided in [11]. We generate an ensemble of series expansion coefficient sets 
randomly distributed according to a multivariate normal distribution, using the fitted central 
values and covariance given in [11]. For each of the sets, we then evaluate the FFs at the three 
q^ values and extract the uncertainties and correlation of these pseudo-data points. 

We then construct a function 


.. ,al) = 


+ X/ ['^LCSR(9fc) “ “m)] (C*] 

ijkl 


djkl \ 
LCSR^ 


-1 


T^J 

-^LCSR 


(Qf) - 


{QI ) ^r, 


ijkl 



( 21 ) 


where are the central values of the pseudo data points of FF i and the corresponding 
covariance matrices (taking into account both the correlation between different FFs and dif¬ 
ferent q^ values). We then sample a likelihood L = e~^ using a Markov Chain Monte Carlo 
(MCMC) approach with flat priors for the series expansion coefficients. From the stationary 
distribution of the MCMC, we extract the central values and covariance of the coefficients. 

Fig. 2 shows the fit result for the B —?• K* FFs in the variable g^. The FF plots, in the 
2 ;(g^)-variable, for the modes B —)• K*, Bg ^ (p and Bg —?■ K* are shown in Figs. 5,6,7 of 
appendix E. The LCSR and lattice pseudo-data points are shown in blue and red. The light 
red dashed band shows the 2-parameter fit from [11]. The solid gray band is our combined 
3-parameter fit result. The numerical fit coefficients, of both fits, are reproduced in Tab. 15 in 
appendix E. As for the pure LCSR fit, the central values, the uncertainties and all correlations 
are provided as ancillary files on the arXiv preprint page (see appendix E). 
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2.0 



Figure 2: Combined LCSR and lattice fit to i? —)■ K* FFs, where lattice data points are indicated in 
red, LCSR points in blue, the gray solid band shows the combined 3-parameter fit and the 
red dashed band the 2-parameter lattice fit from Ref. [11]. In the lower left plot Ti > T 2 for 
< 7 ^ > 0 GeV^ above . 
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We would like to add that the fits are valid under the premise that the LCSR and lattice 
QCD data points and their uncertainties, including correlations, are valid as well as the z- 
expansion being a reasonable model function. There is no evidence against the latter, as we 
have found that adding higher orders in the z expansion and using different parameterisations 
does not change matters. This is an aspect that could change in the future with more precise 
FF determination from LCSR or lattice QCD. Overall the agreement is good. For the central 
values, we find a of 7.0 for B —)• K*, 10.2 for Bg —)■ </>, and 19.5 for Bg —)• K*, for 16 
degrees of freedom. The fits in the z(( 7 ^)-variable, shown in Figs. 5,6,7, are particularly useful 
in judging the quality of the fits. In these figures the FFs times P = {q^ — is plotted. 
The latter serves to cancel the first physical pole, at the resonance R, in the g^-spectrum. 
The remaining slope therefore is a measure of the behaviour of the higher poles or cuts in the 
q'^-spectrum. 

An interesting qualitative feature is the behaviour of the B —)• K* versus Bg ^ (p lattice 
FFs P ■ T 23 {q‘^) (and to some extent also P ■ A 12 ). From Figs. 5,6 is seen that the slopes 
are opposite in direction for the two cases. In a LCSR computation, valid at low such a 
qualitative difference cannot arise since the main difference between the FFs for B —)• K* versus 
Bg ^ <p is due to hadronic input data (which is numerically similar). It is possible that by 
going closer to the hadronic spectrum, at high a more distinct pattern arises in accordance 
with the lattice QCD computation. It will be interesting to see whether this qualitative feature 
which is not yet statistically significant is confirmed in future lattice predictions with higher 
statistics and a more complete treatment of the vector mesons (e.g. physical quark masses). 

To this end we would like to add that differences in normalisation of ||,T (V, Ai and ri, 2 ) 
versus 0-helicity (Aq A 12 and r 23 ) in the LCSR computation are highly sensitive to fy and fy 
decay constants. For instance the 0-helicity FFs depend to ~ 75% on the normalisation of fy 
with the situation being just the opposite for the ||, T-helicity FFs. 

3. Phenomenological applications 

We make use of the updated FFs and their error correlations in predicting experimentally 
accessible observables. More specifically we consider the 5 —)• s flavour-changing neutral current 
transitions (FCNC) B —)• K*, B —)• K*'y, Bg —)> sensitive to physics beyond the 

SM and the branching fractions of the tree-level decays B —)> (p, The latter are of interest 

to extract the CKM matrix element | Wbl and conversely serve as a test of the FF normalisation 
(and shape) when iWbl is taken as an input from other channels and global fits. 

3.1. B —>• at low 

The decay B —)• being one of the golden channels of LHCb, requires no introduction. 

It has received a great deal of attention, particularly in the last decade. By making use of 
the large energy relation [5] , observables have been identified which have reduced uncertainties 
with respect to FFs (e.g [9]). Recent measurements and analyses of several of these observables 
by LHCb [56-59], CMS [60], ATLAS [61] and Belle [62] have revealed a number of potential 
tensions with the SM predictions. Whether or not this is due to new physics or hadronic 
effects is a subject of vital debate [7,8,63-67]. This motivates reinvestigation into predictions 
of hadronic quantities such as the FFs undertaken in this work. 
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In the SM, neglecting the muon mass, the differential decay distribution of i? —)■ K*n 
can be written in terms of of six helicity amplitudes 


hY = N^ + , ( 22 ) 

///^ = iV^CioVW(g2), (23) 


where l = +,—,0 denotes the polarisation of the ii'*-meson. The helicity FFs are 

defined as in appendix D and stands for various corrections to be discussed further below. 
The quantity 


N = HbK: 




3 • 


A1/2 


1/2 


(24) 


is a normalisation factor where Gp stand for Fermi’s constant, e for the electric charge and 
Kb(s) CKM matrix elements. The differential branching ratio is then given by 


dBR{Bq 

dq^ 



E E 

i=±,0X=V,A 


(25) 


Factorisable quark loop contributions are absorbed into Wilson coefficients Cf® and 
which therefore become g^-dependent (e.g. [68] for the definition). The quantities contain 
the NNLL corrections to the matrix elements of the current-current operators [69] as well 
as various “non-factorisable” contributions. The latter entail the effect of weak annihilation, 
the chromomagnetic operator contribution both computed in LCSR [12, 47] as well as hard 
spectator scattering taken from QCD factorisation [70,71]. 

An important contribution arises due to the final state leptons emerging from charm quarks; 
so called charm loops. At high the effect of the broad charmonium resonances measured 
by the LHCb-collaboration [72] has turned out to be substantially more sizeable than antici¬ 
pated [66]. More precisely, for BR(i? —)• R/x/i) the resonance residues are found to be ~ 2.5 
larger with opposite sign from naive factorisation indicating sizeable duality violations [66]. 
It therefore seems well-motivated to discuss the various contributions in some detail. At low 
q^ < 6 GeV^ such effects are thought to be captured in the partonic language of charm quarks 
and gluons. The 0(as) hard vertex corrections [69] factorise in the heavy quark limit into a 
g^-dependent function times the FF. The part which does not contain gluon exchanges between 
the hadron transition and the charm loop factorises non-perturbatively by definition and the 
Q'^-dependent function is given by the vacuum polarisation. The latter can be extracted in 
a model-independent way from e'’"e“ —?■ hadrons-data.^^ These contributions, as mentioned 
above, are included in the central values of the predictions of our work. In addition there is soft 
gluon emission from the charm loop into the R-meson as well as the TC^-meson. Both effects 
have been assessed in LCSR, the former with a R-meson DA [73] and the latter with R^-meson 
DA (for R —7- K*'y only) [74,75]. The combination of the two results is not completely free from 
model-dependence.^^ At q^ approaching the charmonium resonance region, the contribution 
is predicted to be significantly enhanced, rendering the partonic theory prediction unreliable 

^^Cf. [66] for a recent determination. 

^®The problem is that the two effects are computed in two, slightly, different frameworks. It wonld be best to 
compute, in either of the two frameworks, the radiative corrections which wonld then inclnde both effects as 
well as the 0{as) vertex corrections. This conld be a rather challenging as it would seem to require analytic 
results in order to be verify the dispersion relation. 


17 









above about 6 GeV^. These two effects, the soft gluon emission and the charmonium effect, 
can be captured in the region G [0, 6 GeV^] by a linear parameterisation 


Af 






6GeV^ 


(26) 


where a, b are positive numbers and (j)afi are strong phases whose parameter ranges we discuss 
further below. Note that (26) and the replacement of 


Cf 


Cf 


1 + 


q 


6GeV" 


(27) 


in (22) are equivalent to each other. The parameterisation (26) is convenient for low since it 
incorporates the helicity hierarchy A_|_ <C A_ through the FF parameterisation. This results 
in a+ ~ a_ and 6+ ~ 6_. We find^^ 

a±G [0,0.05], 6±G [0,0.2], 

ao G [0,0.2] , bo G [0,0.05] . (28) 


where is mainly fixed at low q'^ by the soft gluon emission [73-75] and 5^ is then determined 
to cover the J/T uncertainty. We vary the phase of the J/T-residue in the dispersion rep¬ 
resentation in the full range motivated by the findings in [66]. Note that the absolute value 
of the residues are known from the polarisation specific branching fractions B —)■ J/^K*. 
The asymmetry between the parameter values of ao,bo and a±,b± is due to the ± directions 
being sensitive to the photon pole (contrary to the 0-helicity direction). At intermediate q^ 
this hierarchy disappears which can for example be seen from the polarisation fractions of 
the B —)• J/'^K* amplitudes or the general result that the helicity amplitudes are degenerate 
at the kinematic endpoint [77]. In summary the uncertainty due to soft gluon emission and 
nearby resonances is covered by the parameterisation (27) with parameter ranges as given in 
(28) and varying the phases (jiafi in the full range. 

Numerical predictions in different q'^ bins for B^ —)• observables (see e.g. Refs. [68, 

78,79] for definitions of the angular observables) and the B^ —)• K *^branching ratio are 
given in Tabs. 6, 7 and 8 respectively. Grucially uncertainties are split into parametric^®, FF 
and non-factorisable charm uncertainties as parameterised in Eq. (26). It is observed that the 
dominant uncertainty of the branching fraction is due to FFs and amounts to about 20% relative 
to the central value. In the case of the angular observables the error is considerably reduced 
by the inclusion of the correlations. Comparing the angular observables ^4^5 with the related 
observables it is noted that the FF uncertainties are comparable. This improvement for 
54^5 observables is due to the inclusion of correlated uncertainties in the FFs. The error due to 
the A-corrections results in comparable uncertainties in both bases. The advantage of using 
fully correlated errors in explicit computation over general parameterisation can be seen by 
comparing the uncertainties in our work versus those of Ref. [8]. 


refer the reader to the appendix of [47] and [76] for recent theoretical discussions of this topic. 

'^®Compared to the parameterisation used in [67[ the value of foo is considerably reduced. For the observables 
presented in this paper this effect has a negligible influence on the values of the uncertainty. 

^®The parametric uncertainties, with values adopted from the PDG [80], include IKbRsI = (4.01 ± 0.10) • 10“^, 
the scale variation fj, = 4.8 ± 0.8 GeV, the &-quark MS mass mb(mb) = 4.18 ± 0.03 GeV and the pole mass of 
the charm quark rric = 1.4 ± 0.2 GeV. 


18 









Parameter 

Value 

Ref. 

TbO 

1.520(4) ps 

[81] 

tb± 

1.638(4) ps 

[81] 

TBs 

1.604(10) ps 

[81] 

AFs/F^ = 2ys 

0.124(9) 

[81] 

\Vcb\ 

4.221(78) X 10"2 

[80] 

\VtbV,l/V,,\ 

0.980(2) 



Table 5; Numerical inputs for the SM predictions. 


For comparison of the B —)• K* observables to existing experimental measurements of 

3 fb“^ LHCb data and the implications for new physics, we refer the reader to [67], where the 
FF results of this work were used for a global analysis of 5 —>• s transitions. 

Values of important parameters used for all the SM predictions are given in Tab. 5. 

3.2. B K*j 

The precise experimental determination of the branching ratio for B —?■ iF*7 provides a good 
opportunity to compare our results for the FFs to experiment. The branching ratio of B ^ K*^ 
is given by 

BR{Bq ^ iF*7) = rB,487r2 , (29) 

where q = u,d. We have introduced the superscript q in addition to the previous section 
because we give separate predictions for charged and neutral modes. The amplitudes are 
related to the limit of the vector helicity amplitudes of i? —)• at vanishing dilepton 

invariant mass, 

„2 

HliB K*j) = lim —HI’'^{B K*i+r) . (30) 

(j2—>-0 e 

They can be written as 

= y yAo(2mfeCf (r±(0) + i2mbA^(0)) , (31) 

where and Aq = A|q2=o = (m^ — is the Kallen-function for the photon final state and 

T±{0) = A(g2)|g2=o which results in T+(0) = 2Ti(0) and T_(0) = 0. The quantity 

is defined in appendix D. For A±(0), the following contributions are included, 

• Corrections to the matrix elements of current-current operators [69]; 

• Hard scattering contributions computed in QCD factorisation [70,71]; 

• Non-factorisable contributions of the chromomagnetic penguin operator Og computed in 
LCSR [47]; 

• Weak annihilation computed in LCSR [12]. 
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V- 





Observable 

bin 


SM prediction 




[ 0 . 1 , 1 ] 

0.897 

± 

0.035 

± 

0.147 

± 

0.050 


[ 1 , 2 ] 

0.436 

± 

0.017 

± 

0.094 

± 

0.014 


[2,3] 

0.400 

± 

0.015 

± 

0.091 

± 

0.010 


[3,4] 

0.409 

± 

0.016 

± 

0.091 

± 

0.008 

10 ? 

[4,5] 

0.432 

± 

0.016 

± 

0.091 

± 

0.010 

dq^ 

[5,6] 

0.461 

± 

0.018 

± 

0.093 

± 

0.012 


[1.1, 2.5] 

0.420 

± 

0.016 

± 

0.093 

± 

0.013 


[2.5,4] 

0.406 

± 

0.088 

± 

0.087 

± 

0.094 


[4,6] 

0.447 

± 

0.017 

± 

0.092 

± 

0.011 


[ 1 . 1 , 6 ] 

0.426 

± 

0.016 

± 

0.091 

± 

0.009 


[ 0 . 1 , 1 ] 

-0.093 

± 

0.000 

± 

0.012 

± 

0.001 


[ 1 , 2 ] 

-0.140 

± 

0.002 

± 

0.036 

± 

0.010 


[2,3] 

-0.072 

± 

0.002 

± 

0.021 

± 

0.020 


[3,4] 

0.010 

± 

0.002 

± 

0.011 

± 

0.026 


[4,5] 

0.085 

± 

0.002 

± 

0.023 

± 

0.030 

^FB 

[5,6] 

0.152 

± 

0.002 

± 

0.034 

± 

0.031 


[1.1, 2.5] 

- 0.122 

± 

0.002 

± 

0.033 

± 

0.013 


[2.5,4] 

- 0.010 

± 

0.011 

± 

0.011 

± 

0.010 


[4,6] 

0.120 

± 

0.002 

± 

0.029 

± 

0.031 


[ 1 . 1 , 6 ] 

0.014 

± 

0.002 

± 

0.011 

± 

0.025 


[ 0 . 1 , 1 ] 

0.330 

± 

0.004 

± 

0.064 

± 

0.018 


[ 1 , 2 ] 

0.749 

± 

0.004 

± 

0.053 

± 

0.019 


[2,3] 

0.825 

± 

0.001 

± 

0.041 

± 

0.009 


[3,4] 

0.805 

± 

0.000 

± 

0.046 

± 

0.005 

Ft 

[4,5] 

0.757 

± 

0.000 

± 

0.055 

± 

0.011 

^ L 

[5,6] 

0.702 

± 

0.001 

± 

0.062 

± 

0.016 


[1.1, 2.5] 

0.782 

± 

0.003 

± 

0.048 

± 

0.016 


[2.5,4] 

0.812 

± 

0.044 

± 

0.047 

± 

0.046 


[4,6] 

0.728 

± 

0.001 

± 

0.059 

± 

0.013 


[ 1 . 1 , 6 ] 

0.768 

± 

0.001 

± 

0.051 

± 

0.006 


Table 6: Standard model predictions for binned 5° —> K*^jjT observables, where the uncertainties 
are split into parametric uncertainties, FF uncertainties, and our estimate of the uncertainties 
due to missing hadronic effects. 
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Observable 

bin 


SM prediction 




[ 0 . 1 , 1 ] 

0.093 

± 0.000 

± 

0.005 

± 

0.003 


[ 1 , 2 ] 

0.005 

± 0.001 

± 

0.009 

± 

0.010 


[2,3] 

-0.096 

± 0.001 

± 

0.015 

± 

0.013 


[3,4] 

-0.163 

± 0.001 

± 

0.019 

± 

0.013 


[4,5] 

-0.206 

± 0.001 

± 

0.019 

± 

0.011 

04 

[5,6] 

-0.233 

± 0.000 

± 

0.017 

± 

0.009 


[1.1, 2.5] 

-0.027 

± 0.001 

± 

0.010 

± 

0.011 


[2.5,4] 

-0.148 

±0.018 

± 

0.019 

± 

0.018 


[4,6] 

- 0.220 

± 0.001 

± 

0.018 

± 

0.010 


[ 1 . 1 , 6 ] 

-0.145 

± 0.001 

± 

0.016 

± 

0.012 


[ 0 . 1 , 1 ] 

0.254 

± 0.000 

± 

0.009 

± 

0.004 


[ 1 , 2 ] 

0.110 

± 0.004 

± 

0.017 

± 

0.020 


[2,3] 

-0.090 

± 0.004 

± 

0.017 

± 

0.027 


[3,4] 

- 0.222 

± 0.003 

± 

0.024 

± 

0.028 


[4,5] 

-0.306 

± 0.002 

± 

0.025 

± 

0.025 


[5,6] 

-0.360 

± 0.002 

± 

0.022 

± 

0.022 


[1.1, 2.5] 

0.048 

± 0.004 

± 

0.016 

± 

0.023 


[2.5,4] 

-0.192 

±0.023 

± 

0.023 

± 

0.023 


[4,6] 

-0.334 

± 0.002 

± 

0.023 

± 

0.024 


[ 1 . 1 , 6 ] 

-0.185 

± 0.003 

± 

0.019 

± 

0.026 


[ 0 . 1 , 1 ] 

0.240 

± 0.001 

± 

0.006 

± 

0.007 


[ 1 , 2 ] 

0.014 

± 0.003 

± 

0.022 

± 

0.025 


[2,3] 

-0.273 

± 0.005 

± 

0.029 

± 

0.042 


[3,4] 

-0.430 

± 0.003 

± 

0.021 

± 

0.031 


[4,5] 

-0.491 

± 0.001 

± 

0.016 

± 

0.020 


[5,6] 

-0.518 

± 0.001 

± 

0.014 

± 

0.013 


[1.1, 2.5] 

-0.070 

± 0.004 

± 

0.026 

± 

0.032 


[2.5,4] 

-0.398 

± 0.022 

± 

0.022 

± 

0.022 


[4,6] 

-0.504 

± 0.001 

± 

0.015 

± 

0.016 


[ 1 . 1 , 6 ] 

-0.358 

± 0.003 

± 

0.022 

± 

0.029 


[ 0 . 1 , 1 ] 

0.653 

± 0.002 

± 

0.009 

± 

0.012 


[ 1 , 2 ] 

0.280 

± 0.008 

± 

0.031 

± 

0.043 


[2,3] 

-0.254 

± 0.011 

± 

0.044 

± 

0.082 


[3,4] 

-0.585 

± 0.008 

± 

0.035 

± 

0.070 


[4,5] 

-0.732 

± 0.005 

± 

0.029 

± 

0.051 

-^5 

[5,6] 

-0.799 

± 0.003 

± 

0.028 

± 

0.039 


[1.1, 2.5] 

0.126 

± 0.009 

± 

0.038 

± 

0.057 


[2.5,4] 

-0.517 

± 0.043 

± 

0.040 

± 

0.039 


[4,6] 

-0.765 

± 0.004 

± 

0.028 

± 

0.044 


[ 1 . 1 , 6 ] 

-0.459 

± 0.008 

± 

0.034 

± 

0.064 


Table 7: Standard model predictions for binned angular observables, where the un¬ 

certainties are split into parametric uncertainties, FF uncertainties, and our estimate of the 
uncertainties due to missing hadronic effects. 
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B+ - 

-)> K*+p,+pL~ 

Observable 

bin 

SM prediction 


[0.1,1] 

0.923 ± 0.036 ± 0.155 ± 0.052 


[1,2] 

0.474 ± 0.018 ± 0.102 ± 0.015 


[2,3] 

0.438 ± 0.017 ± 0.099 ± 0.010 


[3,4] 

0.448 ± 0.017 ± 0.098 ± 0.009 

107 

[4,5] 

0.472 ± 0.018 ± 0.099 ± 0.011 

dq‘^ 

[5,6] 

0.502 ± 0.019 ± 0.100 ± 0.014 


[1.1, 2.5] 

0.458 ± 0.017 ± 0.101 ± 0.013 


[2.5,4] 

0.445 ± 0.095 ± 0.094 ± 0.102 


[4,6] 

0.487 ± 0.018 ± 0.099 ± 0.012 


[1.1,6] 

0.466 ± 0.018 ± 0.099 ± 0.009 


Table 8 : Standard model predictions for the differential branching ratio of —>■ K*^fi~, where 

the uncertainties are split into parametric uncertainties, FF uncertainties, and our estimate 
of the uncertainties due to missing hadronic effects. 



Theory 

Experiment 

10^ X BR(R0 ^ 

3.39 ±0.14 ±0.70 ±0.28 

4.33 ±0.18 

10^ X BR(R+ ^ K*+'y) 

3.33 ±0.13 ±0.72 ±0.29 

4.21 ±0.15 


Table 9: SM predictions and experimental world averages for the branching ratios of 5° —>■ and 

B+ K*+ 7 . The theory uncertainty is split into parametric, FF, and non-factorisable power 

correction uncertainties. 


The first of these corrections is by far the dominant one, leading to a +60% shift in the 
branching ratios. The three remaining ones contribute to the isospin asymmetry (e.g. [12]) of 
which WA is the one which is most sizeable. 

Our predictions for the branching ratios are listed in Tab. 9 along with the experimental 
world averages and are consistent with the experimental results at around la. We would like 
to emphasise that the B —)• K*^ is a FCNC and that the consistency cannot be taken to be 
one to one with a FF measurement. The B —?• {p,u:)ii' decays, discussed in section 3.5, are 
more favourable in this respect. 

Another cross-check is the branching ratio of the decay B —>■ K*e~^e~ at very low that is 
dominated by the photon pole and that has been measured recently by LHCb [82], 

BR{B a:*± 0.2) X 10“^ (32) 


where the superscript refers to An interesting observable is the ratio of this branching 

ratio to the B —?• K*'y branching ratio, since theoretical uncertainties, factorisable or non- 
factorisable, cancel to a high degree. In the SM, we predict 


_ BR(S0 
■^ee,7 — 


^ 7^*0g+g-)30-1000MeV 
BR(R0 ^ K*0'y) 


= (6.3 ±0.2) X 10 


-3 


(33) 
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where the residual error is dominated by FF uncertainties. Combining experimental errors 
in quadrature, from the LHCb measurement and the world average of —)• we 

obtain 

= (7.2±2.1) X 10-^ (34) 

which is consistent with the prediction, albeit with sizeable uncertainties. Finally, for the 
angular observable Fl, that has been measured recently in —)• K*e'^e~ at low [83], we 
predict 

Fl{B x*g+e-^i5-i058MeY ^ q_203 ± 0.003 ± 0.058 ± 0.017 . (35) 

This is in very good agreement with the experimental value, 

Fl{B = 0.16 ± 0.06 ± 0.03 . (36) 

3.3. Bs —>• at low 

The decay channel Bg —t- (j)F^ ijl~ is proceeds via the same quark level transition as R —?■ 
K *and may serve to compare possible deviations. An important difference between the 
two channels is that the iji-meson decays to K^K~, implying that the decay is not self-tagging 
in contrast to B^ K* , where the flavour of the initial R-meson can be inferred from 
the charge of the Ktt decay products of the K*. As a consequence, among the observables 
discussed for B —)■ K*, Ape and cannot be measured at a hadron collider. 

Other than CP asymmetries, the most interesting observables are then the differential 
branching ratio, Fl, and ^ 4 , in the SM and beyond. For these observables, the three pos¬ 
sible sources of difference between the results for Bg —)• and those for B —)■ K* n~ 

are as follows, 

• the FFs are different; 

• differences induced by spectator effects, e.g. weak annihilation; 

• effects due to the sizeable Bg-Bg lifetime difference. 

The FFs have already been discussed in sections 2.5 and 2.6. The spectator effects turn out 
to be very small in the SM and are not relevant compared to the FF uncertainties. For a 
discussion of effects beyond the SM we refer the reader to the appendix of Ref. [12]. The 
lifetime effects are due to the Bg and Bg lifetime difference of roughly 6 % absent for B^- 
mesons. This leads to a difference between the observables defined in the absence of neutral 
meson oscillations, as used in the case of B^ —)• K*, and time-integrated observables, as 
measured experimentally [84]. This difference has to be taken into account when comparing 
theory predictions to experimental data. Details are discussed in appendix G. 

In Tab. 10, we list our numerical predictions for the differential branching ratio and angular 
observables accessible from an untagged measurement Bg —?■ (f)F^^~. The uncertainties are 
treated in the same way as for B —?■ K* . 

3.4. B versus Bg —>• cf)F^F~ 

The similarity of the B —)• K* F ^and Bg —)• ^F^ channels implies that the uncertainties 
of ratios of these observables should be strongly reduced.Theory predicts Bg —)■ (f)F ^to 

this work we have not performed an error analysis on the ratios themselves. The latter would greatly reduce 
the error and could be undertaken if the experimental central values persist with smaller uncertainties. 
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Bs 








Observable 

bin 


SM prediction 





[ 0 . 1 , 1 ] 

1.067 

± 

0.042 

± 

0.155 

± 

0.058 



[ 1 , 2 ] 

0.497 

± 

0.019 

± 

0.099 

± 

0.017 



[2,3] 

0.450 

± 

0.017 

± 

0.096 

± 

0.011 



[3,4] 

0.459 

± 

0.017 

± 

0.096 

± 

0.009 

10 ^ 

dBR 

[4,5] 

0.484 

± 

0.018 

± 

0.097 

± 

0.011 

dq'^ 

[5,6] 

0.516 

± 

0.019 

± 

0.099 

± 

0.015 



[1.1, 2.5] 

0.476 

± 

0.018 

± 

0.098 

± 

0.014 



[2.5,4] 

0.455 

± 

0.017 

± 

0.096 

± 

0.009 



[4,6] 

0.500 

± 

0.019 

± 

0.098 

± 

0.013 



[ 1 . 1 , 6 ] 

0.479 

± 

0.018 

± 

0.097 

± 

0.010 



[ 0 . 1 , 1 ] 

0.311 

± 

0.004 

± 

0.057 

± 

0.017 



[ 1 , 2 ] 

0.732 

± 

0.003 

± 

0.051 

± 

0.020 



[2,3] 

0.813 

± 

0.001 

± 

0.039 

± 

0.010 



[3,4] 

0.791 

± 

0.001 

± 

0.045 

± 

0.006 

Ft 


[4,5] 

0.739 

± 

0.001 

± 

0.054 

± 

0.011 

^ L 


[5,6] 

0.682 

± 

0.001 

± 

0.060 

± 

0.016 



[1.1, 2.5] 

0.767 

± 

0.003 

± 

0.046 

± 

0.017 



[2.5,4] 

0.799 

± 

0.001 

± 

0.043 

± 

0.006 



[4,6] 

0.710 

± 

0.001 

± 

0.057 

± 

0.014 



[ 1 . 1 , 6 ] 

0.752 

± 

0.001 

± 

0.050 

± 

0.006 



[ 0 . 1 , 1 ] 

0.088 

± 

0.000 

± 

0.005 

± 

0.002 



[ 1 , 2 ] 

0.003 

± 

0.001 

± 

0.009 

± 

0.010 



[2,3] 

-0.099 

± 

0.001 

± 

0.016 

± 

0.013 



[3,4] 

-0.166 

± 

0.001 

± 

0.019 

± 

0.012 



[4,5] 

-0.208 

± 

0.001 

± 

0.018 

± 

0.010 

04 


[5,6] 

-0.234 

± 

0.000 

± 

0.016 

± 

0.008 



[1.1, 2.5] 

-0.029 

± 

0.001 

± 

0.011 

± 

0.011 



[2.5,4] 

-0.151 

± 

0.001 

± 

0.018 

± 

0.013 



[4,6] 

- 0.221 

± 

0.000 

± 

0.017 

± 

0.009 



[ 1 . 1 , 6 ] 

-0.146 

± 

0.001 

± 

0.016 

± 

0.012 


Table 10: Standard model predictions for binned, time-integrated Bg —)■ observables, where the 

uncertainties are split into parametric uncertainties, FF uncertainties, and our estimate of 
the uncertainties due to missing hadronic effects. 
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have a higher transition than B —>■ K *which essentially comes from the decay constants 
(cf. Tab. 1) showing this hierarchy. At low and for final state (i.e. = 0) the 

central values of the LHCb results show the opposite effect. 

First, we recapitulate the prediction for the branching ratio of Bg —>■ (^>7 (see appendix A of 
Ref. [12] for more details) versus B —)■ K*'y. The effect is driven by (0)(0) = 

0.89 ± 0.10^®, resulting from the above mentioned decay constants, leads to 


(,) _ BR(R° ^ iF*° 7 ) 
BR(R, ^ 


0.78 ±0.18, 


(37) 


which is roughly 1.5 standard deviations below the LHCb measurement for this ratio, 1.23 ± 
0.32 [85]. Such a deviation can, of course, not be regarded as statistically significant. 

A similar ratio can also be considered for the decay to leptons. 


RK*4>[qi,q2] 


dBR(R0 ^ 

dBR{Bg ^ <p£+i-)/dq\^^g^] 


(38) 


by considering ratios of the differential branching ratios integrated over specified ranges in q^. 
We show a graphical comparison of our predictions using LCSR, lattice and combinations of 
the two for the ratio R^k* to the results of LHCb [56, 86 ] and CDF [87] at both low and 
high q^ in Fig. 3. Again, the results per se are not statistically significant. On the qualitative 
level it is though interesting that the theoretical and the experimental ratio are below and 
above unity respectively. It is hard to see how the theoretical value can move above one, 
through redetermination of parameters, without uncovering a new physical effect. We stress 
once more that we have not undertaken an analysis with correlated errors for this quantity. 
One could easily expect the theory error to reduce down by a factor of two which would result 
in RK*cf)\[ifi] < 1 within uncertainties. We are looking forward to the 3fb~^ results to reexamine 
this issue. 


3.5. |Vub| from B —>• {p,uj)£iy 

FCNC decays such as B ^ K* and Bg —)• are potentially affected by new physics 

and therefore do not provide an unambiguous environment to test FF predictions. The semi- 
leptonic decays B —)■ {uj,p)£^u based on the charged current b ^ u transition occur at tree- 
level, and are therefore less likely to be affected by new physics and serve to test FF predictions. 
In particular in view of the current discrepancies between the B —)■ K* p,~ and Bg —>■ 

branching fraction measurements the normalisation of the FFs per se has become an issue of 
considerable interest. 

The differential branching ratios of these decays, for mi = 0, are given by 

^ (ivwp ^ |V.-)|^ + |V(°>P). (39) 

dBR(H+ ^ p^£+i2) _ tb+ dBR(H° ^ p-£+u) 

d? “ 2 ^ d^ ’ ^ 

where definitions of for i = ±,—,0 as well as A, the Kallen-function, can be found in 
appendix D, with the adaption mx* iTip- The one for B~^ —?• is analogous to B~^ —?■ 

p^i'^u with obvious replacements. 

'^^The central value of this work, which is a more complete update, is (0)/Tj®''“^'^(0) = 0.91. 
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lattice 


LCSR 

LCSR + lattice 

LHCb 

CDF 

0.0 0.5 1.0 1.5 2.0 

(-RjC (#)[1,6] 


LCSR + lattice 


LHCb 


I—.—I -CDF 

0.0 0.5 1.0 1.5 2.0 

{Rk‘ a>)[16,19] 


Figure 3: Our predictions for at low and high using LCSR, Lattice and a combination of 

the two, compared to experimental results from LHCb with and integrated luminosity of 
1 fb-i [56,86] and CDF [87]. 


The most recent measurements of the branching ratios have been performed for B —?■ pii' 
by BaBar [88] and Belle [89] and for B —> ojlv by BaBar [90,91] and Belle [89] respectively. 
We extract [Wbl from the BaBar and Belle data by minimizing the function that reads in 
both cases 


x'(|K.b|)=E[^exp 

ij 





^t'hdKrbl) 


(41) 


where Bl^^ and Bl^ are the experimental and theoretical central values for the branching ratios 
in one q^-hRi and the sum runs over all bins for the charged and neutral mode. Cth is the 
theoretical covariance matrix that includes in particular the correlated FF uncertainties^®. In 
the case of Belle, we use the data up to = 8 GeV^ or 12 GeV^ and take the full covariance 
matrix provided in Ref. [89]. The BaBar dataset consists of a single bin in the low-cp' region 
from 0 to 8 GeV^ and the correlation between the charged and neutral decay is not provided. 
For B —7- piv we obtain 


1 F K 1 — 

1 Belle, g2<8GeV2 

(3.36 ±0.17 ±0.34) x 

lo■^ 

(42) 

1 F K 1 ^^P^^ — 

' “°lBelle, g2<i2GeV2 

(3.25 ±0.14 ±0.34) x 

lo■^ 

(43) 

1 F u 1 ^^P^^ — 

' BaBar, q^<SGeW'^ 

(2.52 ±0.42 ±0.56) x 

lo■^ 

(44) 

and those from B'^ —)• v we get 




iTy 

1'^ublBelle, g2<7GeV2 “ 

(2.49 ±0.34 ±0.32) x 

lo■^ 

(45) 


^®In Ref. [92] , the importance of uncertainty correlations to extract Vub from B —> p£u decays has been empha¬ 
sized. 
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B^plv LCSR + Belle, q^< 8 GeV^ 
B^plv LCSR + Belle, q^< 12 GeV^ 
B^plv LCSR + BaBar, q‘^< 8 GeV^ 
B^UJlv LCSR + Belle, q‘^< 7 GeV^ 
B-^a}\v LCSR + BaBar, q'^< 8 GeV^ 
B^OJlv LCSR + BaBar, q^< 12 GeV^ 
B^ttI V, global fit 
B-^XJ V inclusive, global fit 
CKMfitter indirect 
UTfit indirect 


2.0 2.5 3.0 3.5 4.0 4.5 


I Kb IX10^ 


Figure 4: Our predictions for |Vub| from B —>• piv and B —>• uj^v (blue) compared to global fits to Kb 
from exclusive [80] and inclusive channels [89] and indirect determinations from fits of the 
unitarity triangle [93,94]. 


= (3.25 ± 0.36 ± 0,53) x lO'^ , (46) 

lb.bli.l"'r,Vl 2 G=v== (3-25 ± 0.29 ± 0.46) x lO-® , (47) 

where the first error is experimental and the second theoretical. For the FF B ^ p and B ^ oj 
we have taken into account that it is a 6 —)• n and not a 6 —?■ d transition by scaling the FFs as 
in (C.16). 

Our results can be compared to the value extracted from B —)■ irli/ decays, obtained in 
Ref. [89] from a global fit of BaBar and Belle data to lattice and LCSR computations, 

= (3.41 ± 0.22) X 10"^ , (48) 

or the average of the inclusive semi-leptonic 6 —>• n determinations [80] 

= (4.41 ± 15l(^) X 10"^ , (49) 

where the first error is experimental and the second error is theoretical. Finally we also compare 
our results to the values obtained indirectly from global fits of the CKM matrix [93,94], 

I Kb I CKMfitter = (3-44lQ;gg) X 10“^ , | Kb I UTfit = (3.61 ± 0.12) X 10“^ . (50) 

The various values for [KbI quoted in this section are summarised graphically in Fig. 4. 

The B —>■ {p,uj) FFs do not, and should not, incorporate an S'-wave contribution since the 
{p,oj) —>■ TTvr is necessarily in a P-wave (cf. section 2.3). Hence the experimental branching 
ratios might be too large which in turn leads to a systematic upward shift of j Kb | as extracted 
from these analyses. In Ref. [42] (cf. Fig. 9 of that reference) this effect has been analysed and it 
has been found that the integrated line-shapes of the S'-wave over the interval [mp—Tp, rUp-l-Fp] 
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is around 12% of the corresponding P-wave contribution. This means that if the S-wave is 
neglected altogether then we can expect an upward shift of ~ 6% of the |Vub| values. In the 
BaBar and Belle analysis the S'-wave has not been subtracted systematically. Hence if the 
precision below the 10%-level is to be reached then the experimental analyses have to perform 
an angular analysis^*^ in order to subtract the S-wave. 

Leaving aside the inclusive determinations we conclude that our | Tub |-values from the com¬ 
bined Belle and BaBar analysis are somewhat lower but surely consistent with B —)• ttIv 
determinations as well as the global CKM fits. The values of iVJibl which are considerably 
lower than the average come with large experimental uncertainties and are consistent at the 
level of one standard deviation. The uncertainty is rather large and an updated analysis with 
the full BaBar data set will be more telling. 

4. Conclusions 

In this paper we present an update of the light-cone sum rules FFs relevant for the B ^ V 
transitions Bj^ g —>■ K*, p, u, (j) using new hadronic inputs such as decay constants (appendix C), 
the newly determined light-cone DAs (twist-5) (appendix B) with explicit results given 
in appendix A.4 as well as in form of a Mathematica notebook. To corroborate the robustness 
of our predictions, we have discussed in detail the role of the equations of motion in reducing 
the uncertainties of tensor-to-vector FF ratios and mass scheme dependence in section 2.2 and 
appendix A.5 respectively. The impact of the H-meson being an unstable particle is analysed 
in section 2.3. 

An important point of this work are the easy-to-use numerical expressions of the FFs, pro¬ 
vided to the phenomenological community, that allow to retain all the uncertainty correlations 
among the FFs in phenomenological analyses. This is of particular importance for predicting 
angular observables that involve ratios of FFs. A ^-expansion fit, Eq. (16), to the FFs is 
provided along central values, uncertainties, and correlation matrices for the expansion coeffi¬ 
cients; available in electronic form as ancillary files on the arXiv webages (see appendix F for 
details and Tab. 14 for the central values). The parameterisation is chosen to transparently 
fulfil the two exact relations among the FFs at = 0. In addition we performed combined 
fits to LCSR and lattice computations of the FFs. This serves on the one hand to obtain 
predictions for the FFs valid in the full kinematic range, on the other hand as a cross-check 
of the consistency between the two complementary approaches, as they have to coincide for 
intermediate values; good agreement is observed between the two. Likewise the z-expansion 
coefficients and the correlated uncertainties of the combined fits are downloadable as ancillary 
files (and central values in Tab. 15). 

A phenomenological analysis is performed using the updated predictions and a new treat¬ 
ment of theoretical uncertainties. In sections 3.1 and 3.2, we have given updated Standard 
Model predictions for B —)• K*pL~ and B —)■ iF *7 observables, taking into account LCSR 
calculations of several hadronic effects beyond FFs and an estimate of the uncertainties due 
to missing hadronic effects, notably contributions from charm quarks. Potentially relevant 
long-distance effects which have not been fully computed are the complete 0{as) charm loops 
effects in one single framework as well as weak annihilation at 0{as)- The latter could be 
sizeable since the chiral suppression might be relieved for radiative corrections. Likewise, in 

^°For the tttt final state one can make use of the isospin in order to dednce whether it is an S- or P-wave 
contribution [38]. 
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section 3.3 we have given predictions for Bg —)• observables, showing (in appendix G) 

that the effect of the finite Ss-lifetime difference is negligible in the Standard Model. Our pre¬ 
dictions are particularly relevant in view of several apparent tensions between Standard Model 
expectations and experimental measurements observed recently. A crucial question to address 
in the near future will be whether these tensions are due to underestimated hadronic effects 
or physics beyond the Standard Model. Our improved FF predictions can play an important 
role in answering this question. 

In section 3.5, the new FF predictions were used to extract the CKM element |Fub| from 
BaBar and Belle measurements of B —)• {p,uj)iu decays. Barring some disagreement among 
the experiments, we find good agreement of our predictions with other exclusive predictions, 
e.g. B —)• and global fits. Our results contribute to the enhancement of the ongoing 

tension between the exclusive and inclusive determination of |Vub|- From another viewpoint 
the encouraging agreement with other exclusive channels serve as a test of the FF normalisation 
which might become particularly important if the disagreement of B —)> versus Bg —)• 

becomes more significant. Our predictions for |Vub| have a relative theory uncertainty 
at the level of 10% from B ^ p and 12-15% from B ^ oj, showing the potential of future, 
more precise measurements of these semi-leptonic decays to improve the precision on |F(ib|- 

We conclude by emphasising that our improved FF predictions are important for the tensions 
in both 6 — 7 - s channels and the determination of | Wb | • These questions can be further examined 
with future experimental data to which we look forward. 
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Addendum 

Significant changes from vl to v2 of this paper are summarised here. In appendix C the 
extraction of the (longitudinal) decay constants fy has been redone, following earlier analyses 
in [47,74], using updated experimental results and more refined definitions. This led to changes 
in the decay constants, and therefore FFs, of around 2-3% depending on the vector meson. The 
twist-4 DA A||, differing from previous results in the literature, and the twist-5 DAs DAs 
have been determined and implemented into the calculation cf appendix B.l. Using the latter 
the tree-level contribution to the FF have been recomputed. This led to downward changes 
of around 4%, 6 % and 10% for the ±-helicity, Ao,i 2 and T 23 FFs respectively. The EOM in 
(6)-(9) have been explicitly checked at the tree-level including in particular the new twist- 
4,5 contributions. Explicit results are given in appendix A.4 and in form of a Mathematica 
notebookBSZ .nb. The effect of the finite Bg-Bg lifetime difference discussed in appendix G is 
now included in all predictions of Bg —)• cjyp^p~ observables. 
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There are only minor changes from v2 to v3. In particular the numerics remain unchanged. 
Changes include clarifying remarks in the text and the adaption of the notation of the twist-4 
DAs in accordance with the literature (cf. footnote 29). 


A. Aspects of the LCSR determination of the form factors 

A.l. Equation of motion and correlation functions 

The LCSR FFs are computed from a correlation function of the type 

C[F] = i j {K* {p)\T{sTh{Q)JB{x)]\^) , (A.l) 

where F is the Dirac-structure and Jb = rnbbi'y^q an interpolating field for the R-meson. In 
fact, the projection on the R-meson through a dispersion relation and the Borel transformation 
can be seen as the substitute of the LSZ-formalism for the R-meson. It is well-known that at 
the level of correlation functions EOM (5) are corrected by contact terms. This results in 

q‘'C[icr^^{'y 5 )] + - (m^ ± rnb)C['y^{'y 5 )] + C[i{d + d)/,( 75 )] - 2 iC'[D^ ( 75 )] , (A. 2 ) 

(5) 

where ^ denote the contact terms. A heuristic derivation of the contact term follows from 
the time derivative acting on the time ordering of the operators which leads to a commutator 
expression 

Ajf) = - J {K* {p,r])\[siaf,o{-f5)b{0), Jb{x,0)]\0) . (A.3) 

Using the canonical equal-time commutation relation for the 6 -quarks, { 6 L(x, 0), 6 ^ 3 ( 0 )} = 
d^^\x)6ap, leads to^^ 

A^ = 0 and A^ = -imbfyr]*^ . (A.4) 

The crucial point is that the contact term is a local term which does not affect the extraction 
of the FFs at all since it does not enter the dispersion relation. Hence the FFs which are 
determined from the correlation function in LCSR obey the EOM. More precisely the EOM 
impose constraints/correlations on Borel parameters and continuum thresholds of sum rule 
parameters. 


A.2. Correlation of continuum thresholds and Borel parameters 

Each correlation function obeys a dispersion relation. Using the notation = Ci [ 7 ^]Rf 

Cl [ 7 ^] can be written as 


C'i[7^] 



ds 


pib^Ks) 
s — pg — iO 


m%Vi{qbfB ^ 

p%-m% + iO Js^ s-p%-i0 ’ 


(A.5) 


where Sc marks the continuum threshold and “cut” stands for the beginning of the discontinuity 
which is just below = m^. Since (A. 2 ) is valid for any p^ it follows from the representation 

the computation we have assumed that the vector meson is at rest. The result (A.4) is the covariantised 
version. Alternatively we could have derived the contact term directly from the path integral through field 
transformations or the (covariant) r*-product. Both of which should directly yield a covariant result. 
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(A.5) that the EOM are valid for the densities pi{s) point by point, i.e. locally. In particular 
subtracting the FF FOM (6) from (A.2) we obtain 


Js, S-p%-l^ 


-{rus + rub) 



ds 


Pi[7^](s) 

s — pg — zO 



Pi\Dp]{s) 

s — — iO 


(A.6) 


for the direction with somewhat elaborate notation. This is of course true for the exact 
density as well as for the density ^lc-OPE computed from the light-cone OFF. The semi-global 
quark-hadron duality, or sum rule approximation, consists of replacing the integral on the 
right-hand side of (A.5) by: 



ds- 


Piis) 


s — p^ — iO 





pLC-OPE(g 

- p| - zO ’ 


(A.7) 


where Sq is some effective threshold parameter which is expected to lie somewhere between 
(niB + 2m„) ~ 30.9 GeV^ and (ms + zzzp)^ ~ 36.6GeV^. A simple way to achieve consistency 
with the FOM (A.2) is to impose . From a physical perspective this is a natural 

choice since the currents are of identical quantum numbers and therefore couple to the same 
spectrum of states. Below, it is argued that the FOM strengthen this point implying a high 
degree of correlation of the continuum thresholds. 

Our main point is that since Pi <C Ti,V, which we infer from the closeness of (11) to 
unity (cf. Fig. 1), a relative difference between and can only be compensated by a much 
larger change in The latter corresponds to a gross violation of semi-global quark hadron 
duality which we exclude; partly on grounds of past experience with LGSR. 

Let us illustrate this more quantitatively by considering r-ratio which are not accidentally 
close to zero cf. Fig. 1. Let us choose r_L(0) ~ 0.957, for = 35 GeV^ and ±0.5 GeV^, 

with fixed Borel parameters, the EOM (6) requires = (35^®) GeV^ which are considerable 
shifts. This corresponds to a change in the FF Pi of roughly 50% in both directions. The 
situation is similar for the other directions.From this we infer that a difference of 1 GeV^ 
on the two continuum thresholds Sq^ — is at the upper boundary of what seems plausible. 
For = 35(2) GeV^ this can be imposed by correlating the two FFs by 7/8 (i.e. 87.5%). 
The same line of reasoning applies to ry and rg+t (H)- Yet for ro+t the numerics are less 
compelling (cf footnote above) and we restrict the correlation between to 50%. There are 
two further correlations at = 0, namely Ti(0) = r2(0) which is of the algebraic type and 
Ao(0) = ^ 3 ( 0 ) which is required to avoid an unphysical pole at q^ = 0. This leads to Sq^ = 
and Sq° = Strictly speaking the latter two are only exact at = 0 but since we refrain 
from assigning a g^-dependence to so the relation is assumed throughout. 

In summary the following correlations are imposed, 

corr(so^s^) = 7/8 , corr(so^ Sg i) = 7/8 , corr(so Sq^^) ^ 1/2 , (A.8) 


^^With respect to vl of our paper the argument is even stronger as the twist-4 contributions do all satisfy the 
eom. 

^®The main conclusions remain unchanged when other points are chosen. For example r||(10GeV^) ~ 1.022 
requires shifts of = 35 (+ 10 ®) GeV^ for Sq^ = 35 ± 0.5 GeV^ again corresponding to shift of around 50% 
of I?2(10GeV^). In fact the -1-10 GeV^ is only compatible within 10% with the EOM. The EOM are only 
satisfied asymptotically for this case. The situation for ro+t is less compelling and requires shift of ±2 GeV^ 
for an original deviation of 0.5 GeV^ which is still remarkable though. 
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and the full correlations corr(s^\ = 1 and corr(s^°, = 1 together with (A.8) imply 

corr(so\so^) = 7/8 , corr(so^ sj") = 7/8 , corr(so Sg °) = 1/2 . (A.9) 

The reader is reminded that we have argued in section 2.4 for a correlation of the type 
corr(s^, Sg®) = 1/2 where F stands for any FF and Sg® is the continuum threshold for the fs 
sum rule. 

So far we have not discussed the role of the Borel parameter. In principle one could argue 
that the Borel parameter and the continuum threshold could conspire to satisfy the EOM. 
Whereas this does not seem to be very viable from the point of view of physics it is in addition 
not credible on grounds of the actual numerics. For example doubling the Borel parameter 
of the light-cone sum rule keeping the Borel parameter of the Jb sum rule fixed, 

leads to a change in the FFs Ti and V of just one percent. Doubling the sum rule parameter is 
outside the validity range since it enhances the continuum contributions relative to the B-pole 
contribution. For example for ri(0) the continuum contribution becomes 42% by doubling 
Hence the Borel parameter cannot balance a change in the continuum threshold of sq 
of 1 GeV^. Hence it is legitimate not to enter the in the discussion. The sensitivity 

of the /b sum rule to the Borel parameter is slightly higher presumably because the 
local condensates are more vulnerable to (semi-global) quark hadron duality violations. This 
uncertainty is important for the FF prediction per se but only enters the EOM by a global factor 
and is therefore not relevant for the discussion of this section. Hence, we fully correlate the 
uncertainties of the Borel parameters corr(M|^, Mp) = 1 which is justified since the variation 
in So are responsible for the bulk part of the uncertainty. It should be added that it is the 
variation of the parameters and sq that addresses the validity of the semi-global quark 
hadron duality. 

A.3. Remarks on the explicit verification of the EOM at tree level 

In view of the importance of the EOM for the determination of the ratio of tensor-to-vector 
EFs we discuss in some more detail the DAs that enter the EOM and the consistent handling 
of the projection of the correlation functions on the EFs. We have explicitly verified the EOM 
at the tree level for all five structures appearing in (A.2) including the derivative FFs as well 
as the strange quark mass terms at twist-3. The strange quark mass terms cancel non-trivially 
between the explicit term in (5) and ©(m^j-correction in the DAs. More detail on the latter 
can be found in section B. The mechanism that guarantees this interplay of light and heavy 
quark mass are the EOM of the light DAs. Explicit results are given in the next subsection A.4 
and in form of a Mathematica notebook. In addition we have verified that the renormalisation 
of the composite operators (cf. section A.5.1) are compatible with the EOM as expected from 
first principles. We remind the reader that the completely tractable issue of contact terms has 
been discussed and resolved in section A.l. Let us add that the covariant derivative between 
the strange and beauty quark can be exchanged by using the following algebraic identity 

C\i{d + d)fi] — 2iC\Pij] = —C\i{d + d)fi] + 2iC'[D^] . (A.10) 

To this end we would like to discuss the consistent handling of the projection onto the 
structures (2). First, we note that for the EOM to be satisfied the projection on the 

Lorentz structures ought to be handled consistently for all structures appearing in (A.2). Eor 
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P ^3 extra care is in order, see e.g. discussion in [47], since q-{p+PB) = p's — n^v equals m\ — my 
only if p^Q is on-shell. Hence in the computation the following projectors ought to be used 

P2 = i{<l-{P + PB)ri*^-{v*-q){p + PBY] , 

p'^ = i{'q*-q){q^ -- Ap + PbY} , (A.ll) 

q-ip + PB) 

which we denote by a lower case p. The important point is that these projectors are transverse 
q ■ Pi = 0 even in the case where (off-shell). We should add that the actual effect 

on the standard tensor and vector FFs due to the difference of using either P 23 ^23 

rather small (numerically around 2 %) since the the sum rule aims at imposing ~ by 
construction. The latter might be taken as a measure of the quality of the sum rule. See also 
the discussion on the optimisation of the Borel parameter given in appendix A. 6 . 


A.4. Explicit tree level results 

The FF densities Rt^ ^ g^Vi 230 given in appendix A.4.2 relate to the FFs as follows 


F{qY 


cp - ^ - B 

IBbIb 


[Rf{u, qY) du 


(A.12) 


with 

Y 

CTi 2 3 = i- ) ^v,Ai — ~{b^b ± mv) , CAi 0 = =t- ) (A. 13 ) 

’ ’ ’ my 

being a matching factor, with a slight abuse of notation, to translate from {V, Ai, A^, Aq} -h- 
{Vi, V 2 , V 3 , Vp}, cf. (3). The symbol B denotes the subtracted Borel transformation explained 
in the next subsection. As an example we write 


V{qY = -(mB + my) 


mb 


m%fB 


B 




du 


fy 


= {mB +my)^^^ 


du 


LJO 

2 


2A 

v2'i\ / a ^2 


(t)±{u) + 0{my) + ... 


mifB 


'UQ 


))/M + O(m^) + ... , (A.14) 


where the dots stand for contributions from other DAs and h{u, qY and uq are defined in 
Eqs. (A.17) and (A.18) respectively. Simpler and more definite expression for the DAs (p±,\\j 
can be found in vl of this paper which fully agree with the current version. Those results 
agree with the expressions given in [ 68 ] for Ti( 0 ), V{qY and ^ 0 ( 5 ^) but differs slightly for 
Ai(p^) due to the previously discussed handling of the projections (cf. section A.3). In practice 
the results for each FF are numerically small but for this work is of importance since we are 
interested in a precise determination of the ratio of tensor-to-vector FF. For example Ai differs 
byaprefactor {m\—my)/{q-{p+pB)) = {fn\ — m‘y)/{p\ — m'y) = {m‘^ — uqY/{umY) + 0{m‘y). 


A.4.1. The subtracted Borel transformation 

The correlation functions of the explicit tree-level results are given by integrals of the type 


In = / du 


Bju) 


A = ml — up% — uq^ 


(A.15) 
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The subtracted Borel transformation 13 is defined as the two-fold operation of taking the Borel 
transformation B followed by the so-called continunm subtraction. From B[Ii] the higher 
B[In>i] can be found by 

Hence, if B[Ii] is fonnd explicitly then the problem is solved. To do so one can write 


h = 


^ A 


du F{u) 2 n 1/2 - 2 n 

/ -T’ Hu,q ) =-{rrih-uq ) , 

/o u h[u) -Pq u 


(A.17) 


and apply the standard Borel transformation B[l/{m‘^ ~ Pb)] ~ ® subtraction 

(corresponding to cntting of the dispersion integral at s = sq) is implemented by imposing 


u> Uq 


B[Ii] 



h(u,q‘^) 

F{u)e 

- du , 

u 


Uq = 


ml — q^ 

So - 


(A.18) 


Using formula (A. 16), taking into account the rub-dependence of uq and h, results in 


B[h] = 


h(u,q ) , _ ®0 

F(u)e F(uo)e 


luo 


M^u'^ 


Uq (so - q‘^) 


B[h] = X 


/ 1 / ^ H{^u,q 

1 / /’^ F{u)e 


luo 




du + 


_ 'iQ 
e aP” 


ul (so - q^) ^ 


(F(uo) (1 -k Xso) - UqF'{uq)) 


BlUl = ^ 


h(u,q'^) 

^ F{u)e Af 2 


UO 


M^u^ 


du -k 


_ SQ 

e aP 


^^0 (so - 


(A.19) 


{F{uq){2 + 2xso + ^so) - uoB'(uo)(2 -k J -k noF"(t6o)) ) , 


where Xgg = (so — q^)/M‘^. 

To this end we wish to comment on the different techniques of Borel transformation. Method 
i): substitute explicit DAs and then integrate over the DA-parameters to obtain analytic func¬ 
tions in (and q^), take the discontinuity, obtain the dispersion relation and then perform 
the Borel transform on the dispersion relation which corresponds to the standard Borel trans¬ 
formation. This method has been pursued for example in [2] for the radiative corrections. 
Method ii): using partial integration rewrite the integrals over the DA-parameters such that 
they take the form of a dispersion relation and then apply the standard Borel transformation. 
This method has been applied in [68] and the vl of this paper to present tree-level results of 
a few DAs. Method iii): the one described above and used here. Methods ii) and hi) have 
an advantage in that one can substitute other DAs after performing the Borel transformation. 
It should be added though that with method ii) that for performing the Borel transformation 
assnmptions on the endpoint behaviour of the DAs might have been made. 
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A.4.2. Explicit tree-level correlation functions 


The relation of the correlation functions Rp (cf. also Mathematica notebook notebookBSZ .nb) 
given below to the FFs is given in (A. 12) and A = — uq^. The definition of the DAs 

and the actual form and values chosen are described in section B. We find (with Rpp = 0 by 
definition) 


rrihfy {uurriy + A) (j)±_ umi,myfyh^^^^ ^ ^ 




i?Ti = 


2 A 2 


2 A 2 


2A 


/f my (mf + A + q^) ^ uflmygy 


8 A 2 


4A3 


A 2 

^( 1 ) 


2A 


. 2 \ ^( 1 ) 


Rt. = - 


umbruyfy [q^u - A{u - 2)) + q^uml + m^) ^ fy'mv {ml - q^) dj 


2 A2 {q^u + m2) 2 2A {q^u + ml) 

{q^ {q^u^ + Aq^uu + A^u) + q^ml {Au — q^ {u^ — 4m + 3)) + q^(3 — 2u)ml — mf) mfmyfyAp 

4A3 {q‘^u + ml) 3 

/fmy {ml - q^) ^ {q^u + m^ + A) ufymy {ml + q^) g^q 


+ 


Rt. = - 


8A2 {q^u + ml) 2 2A {q^u + ml) 

mhmyfy {q^{u — 2)m^ + m^ + q^ {q^{—u) + q^ — Am)) 

A 2 {q^u + ml) 2 

mhfy {q^{u — 2)ml {uumy + A) + m^ {uumy + A) — q'^ {q^u {uumy + A) + Au^my)) dj_ 

2A2 {q^u + ml) 2 

{q^ {q^u - A(m - 2)) + q^uml + mf) umbmyfyh)^^ ^ {ml - q‘^) /|imy<hgl^^ 

“I” ^ A 7 oZ ^ ^ “1“ 


2^2 ^ 2 ^ 2 
m-b - (f) ^ {(fu + ml + A) fymyg- 


2A {q'^u + ml) 


+ mlmyfyAp X 


+ 


8 A 2 {q^u + m 2 ) 2 

{q^ {q^u^ + Aq^uu + A2m) + g^m^ (Am — q^ (m^ — 4m + 3)) + g^(3 — 2u)ml — mf) 

4A3 {q‘^u + ml) 3 

{ml + g^) ufymygl^ _ {q^{u - 2)m^ + m^ + g^ {g^{-u) + g^ - Am)) mbm^ 

2A {q^u + ml) A^ {q^u + ml) 2 

{q^{u — 2)ml {uumy + A) + m^ {uumy + A) — g^ {q^u {uumy + A) + AM^m^)) mbfy4>± 

2A2 (g^u + m2) 2 


+ 


+ 


fy {uumy + A) (/)_L mbflmygl m^ {2ml + fv^-i- 

2 A 2 " 4A2' 8A3 

uml^fl? {ml {2q^u + 3A) + Ag^M + 2m^) umbfymygl 

2A2 {q^u + ml) 2 A {q^u + ml) 
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rriyfy {q^{u — 2)m^ + [q^{—u) + q^ — Att)) umyfyh\^^'^ 

i^q^u + m^) 2 2A {(pu + ml) 

[ml {q‘^{u — 2){uumy + A) + 2Au^my) + m^ {uumy + A) — q^{q^u {uumy + A) — 

A{u - 2)u^ml)) ^^^ {p^ + ml)^ ^ ^ + 2w - 3) + 3A\) - 

Aq^u {(pu + Au) + q^ml ( 2 ( 7 ^ (it^ — 4it + 3) — 3A) + mf [A + g^(4tt — 6)) + 2m®) x 
ml.fpA± 

8A3 [q'^u + 777,2^ 3 ’ 

[q'^u + mf - A) 3mf/ym^^ll^^ + mf) myfy^ 

2A*q^u 2AV 

{—ml {2uumy + A) + (—n) {2uumy + A) + A (A + u{2u — l)m‘y)) mhfymv^^v'* 

A^q^u 

{-2q^u - 2ml + A) mhfvmygf^ ^ {uumy + A) mbfymvg± 

AV + ^2^2 + 

{2q'^u + 2q‘^uml + 2mj — 2A^ + Aq^{2>u — 4)) m^/y 

A^q^u 

{2Aumy — q^ {uumy + A)) fy<j)± {ml + 2A + q^{—u) + q^) myPp^fy 

2AV 2AV + 

mf, (2mg + A) fpA± mf/|',mf.G)[ 

8A3 2A3g2 


mbfymv {ml {2uumy + A) + q"^ (— {2uumy + A)) — A (A + u{2u — l)my)) ^y^^ 

A^q^u 

mbfymv {uumy + A) gp 3m3/ymy {—ml + A + yl|| ^ 

A 2(72 2A^g277 

"^y^±,3/y i^b + 2A - m6/||.m^ { 2 ml “ A + 2g2) 5,(2^ 

2AV A| 33^2 

2myfy {—2q^ml + m^ — A^ + + q^) fv^\\^ 

A^q^u 2A‘^q^ 

mlflmfyG\ _ uml,fp4>± 

2A3g2 Ag2 ’ 


-R©! = 


4A 
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Rt>^ 


Rt^p 


flmv {q^ - ml) ^ 


A [q'^u + ml) A [q‘^u + ml) ’ 

(m^ (^6q^u + 2Aq‘^u + A^) — A (A^ + Ag^(5 — 3u) — 4q^u) + 6mf (A — g^(n — 2)) — 6mf) 

4:A^q‘^u 


3 ( 1 ) (m^ {uumy + A) + q^u {uum'y + A) — Au’^m'y) fymyg^ 




A 2 g 2 


+ 


flmlr {ml (V^ + A) + SAq'^u + 4m^) ^ m^ml,/^ {q^u + ml) /ifj 


A2g2 


+ “ 2) {2uumy + A) — 


{2q^uml — A (a — 2q^u) + 2mf) fymyGj 

4AV 

A {uumy + A)) + mf {2uumy + A) — uq^ {2uumy + A) + A^u^my + 
Aq^ (A + u{2u - l)ml - Au) ) , 


+ 


{ml {uumy + a) + (— {uumy + A)) — Au^my) fymvg\ 

A2g2 

{—ml {2q^ {2uumy + A) + A {uumy + A)) + mf {2uumy + A) + q'^ {2uumy + A) + 
Aq‘^ (A + u{5u - l)my) + A^u^my) ^ 

{—Qmf (A + 2g2) + 77^2 (^_a2 + Qq^ + lOAg^) + 6mf + A (A2 — 4q^ — 3Ag^)) fymyA^^'^ 

AA'^q'^u 

{ml + 2A- g2) _ (Amg + Amf + (5A - V)) ^ 

A'^mb A3q2 

{ml — q^) my/y/ly 4 {—2q^ml + ml — A^ + q^) my/y 
A'^q^mb A^umb 

{-2q^ml + 2ml + A (2g^ - A)) _ 2nm^/^0_L 

4A3g2 Amb 


A.5. Scheme dependence of the form factors 

In many determinations ol B ^ V,P FF calculation in LCSR the pole mass scheme is assumed 
to be the appropriate scheme for the 6-quark mass. For i? —?• vr FFs it has been found that a 
conversion to the MS-scheme leads to minor changes only [20,53]. The explicit appearance of 
mb in the EOM (6)-(9) deserves a reinvestigation of the issue of scheme dependence. 

In LCSR calculations one distinguishes between a factorisation scale g?p ~ — ml ~ 

©(mbAgco) and a renormalisation scale /iuv = mb- The former is the separation scale of the 
LC-OPE and the latter is the scale of the composite operators e.g. the tensor or vector bilinear 
quark currents. For the analysis in this appendix, and throughout the paper, we adopt the 
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B ^ K* 

//2[GeV2] 

ri(o) 

V(0) 

Ti(0)/V(0) 

pole 

4.8 

0.282 

0.341 

0.828 

MS 

4.8 

0.271 

0.330 

0.821 

MS 

8 

0.293 

0.349 

0.840 


Table 11: As mentioned in the text the tensor FFs are understood to be evaluated at the scale /iuv = 
mb by one-loop renormalisation group running. Note — mf ~ 4.8 GeV^ is the 

standard factorisation scale of the LC-OPE used throughout. The values are for central 
values of the input parameters and differ slightly from that obtained from the Markov Chain 
Monte Carlo. 


strategy to lower ^uv to //ir in the actual computation and then use renormalisation group 
running to scale the tensor FFs from to Ti(q‘^)\^^y=rnb- This makes it clear how 

the EOM are obeyed at any step of the computation. More details on the renormalisation of 
the composite operators are given in the next section. 

One can switch back and forth between the pole and MS-scheme by replacing = 

mb{iJ.m)Z^^ = mb(/rm)(l -h ^^0^(4 - 31n(mft///^)) -h 0{al)) (with Cp = 4/3 in 

QCD) in the tree-level computation and expanding to first order in The additional scale 
IJLm is introduced (through the MS -scheme) for the same reasons as /rj^-scale mentioned above. 
In Tab. 11 examples of FF determinations in both schemes are given. We infer that the 
impact of changing from the pole to the MS-scheme for the FFs is around 4% which is sizeable 
but controlled. Yet the ratio of FFs changes by only 1% which is rather small and therefore 
substantiates the robustness of the tensor-to-vector FF ratio which is one of the main points 
of this paper. The //-dependence entering through the //-dependent MS-mass is reflected in 
the pole scheme through a larger uncertainty in the pole mass itself; = 4.8(1) GeV as 

compared to fhb{rhb) = 4.18(3) GeV [80]. 

A.5.1. Renormalisation of composite operators and compatibility with EOM 

The aim of this section is to clarify the renormalisation of the composite operators entering 
the EOM (5) with particular focus on the mb quark mass. The following shorthand notations 
for the operators 


Oi = Od = 2 siDbi b , 02 = Oqt = id''{sia , 

O3 = OmV = {ms + mb)s'y^b , O4 = Oas = idf,{sb) , 


is introduced. The mixing matrix is defined by 


of) = ZgZijOj , 


(A.20) 


(A.21) 


where Zq is the external leg or wavefunction renormalisation. Through an explicit computation 
it is found that 

/2 2 6 6 \ 


Zij — Sij “t" Cp 


(Xa 


1 


Ap 4: — d 


0 0 0 0 
0 0 (2-6) 0 
Vo 0 0 8/ 


(A.22) 
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It is noteworthy that the renormalisation of the operator requires the additional diagrams 
where a gluon originates from the vertex through the covariant derivative. The operators 
OmV,dT,dS do renormalise multiplicatively since they are of lowest dimension (effectively three) 
and differ in quantum numbers when the contraction of the total derivative is undone. The 
operator Od is of dimension four and the dimension three operators can and do mix with Od- 
In the notation (A.20), the operator identity (5) reads 

Od + Omv + OdT - Ogs = 0 . (A.23) 

It is readily verified that the renormalisation (A.22) is compatible with the EOM (A.23). As an 
additional check let us mention that from the diagonal elements Zg-diag(Z) = {Zd, Zqx, Z^Zy, Zqs) 
one infers Zs = Zg^ = 1+6A, Zy = 1 and Zt = Zgx = 1 —2A with A = Zm = 1 —6A. 

From the latter the well-known anomalous dimensions 7 ^*^^ = —BCp, = 0 and 7 -^^ = 2Cf 
of these operators follow (notation: 7 ^ = 7 ^^^ + C*(as))- 

At last we turn to the issue of the impact of the mass renormalisation on the composite 
operators. From the mixing of operators in (A.22) it is clear that the renormalisation of Op is 
affected by a mass scheme change. This can be seen by writing somewhat symbolically Zi^ = 
ZD(mV) = ZraZov- So in Summary going to the pole scheme enforces a finite renormalisation 
of the operator Op since changing from MS to the pole scheme corresponds to a finite shift in 
the ratio of the Z^-factors. Most importantly the renormalisation of the composite operators 
Oy, Ot and Os, on the other hand, is not affected by the mass scheme. Hence it is legitimate 
to use the MS-scheme to renormalise them. This is fortunate since the Wilson coefficients are 
evaluated in the MS-scheme and together this guarantees the cancellation of the /ruv-scale 
between the Wilson coefficients and the matrix elements. The scheme independence of the 
operators Oy, Op underlies or partly explains the small changes in the FFs when going from 
the pole- to the MS-scheme (cf. Tab. 11). 

A.6. Remarks on fixing the Borel parameter 

A sum rule for a FF F{q‘^) of a process B ^ P, F-transition may be written as 

= / pF{s,q^)e~PP~ ds , (A.24) 

Jml 

where is the Borel parameter. The goal of this section is to show that two seemingly 
different methods for fixing are equivalent. For this purpose we introduce the following 
notation 

{x{s))q2^M2^ / x{s)pF{s,q^)e ds . (A.25) 

Jml 

We note that F{q^)M'^ = The two methods are: 

• extremising the Borel parameter: If one were to succeed in computing the sum rule 
exactly, which would imply^^ 

Pf{s, q^) = 6 {s - m\)F{q^) + 0(s - Sc)(Tf{s) , (A.26) 


^For the sake of illustration we employ the narrow width approximation which is justihed for the B-meson. 
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then Eq. (A.24) would remain valid for any Borel parameter. In practice the partonic 
evaluation through the OPE is optimised by using a large Borel parameter, with just the 
opposite being true for the projection on the lowest hadronic state. Hence a compromise 
value has to be found, ideally in a region where shows an extremum in M^. 

This is imposed by 


(A.27) 


daughter sum rule in m\: One may write a daughter sum rule for as follows 


X _ ('S)q2,M2 


(A.28) 


Note that using Eq. (A.26) satisfies (A.28) exactly as it should. 

It is readily seen that Eqs. (A.27,A.28) are the same and hence the two methods are equivalent. 


B. Light-cone distribution amplitudes 


B.l. Distribution amplitudes including Ay and the new twist-5 DAs 

A general review on the subject of LCDAs can be found in [95] which is by now over thirty 
years old. The main concepts for the vector DAs are explained, in some details, in the more 
modern write-up [96]. 

Light-cone physics is conventionally discussed by introducing two light-like vectors say z 
and p (i.e. = p^ = 0).^® The close to light-like separation x and the meson momentum p 

{p^ = my) can be expressed as linear combinations of 2 ; and p 


z^ = Xf,- p^^{xp - J (xp)2 - x^ml,) , Pf,=p^- > 

my v 2 pz 

whereas any vector such as the vector meson polarisation r](p) is decomposed 


(B.l) 




pz vp 

^ Pa i" . 
pz pz 




(B.2) 


into r/-*- and the p and z direction e.g. [96]. Above pz = p ■ z etc is understood. 

The rigorous definition of the LCDAs is given for quark bilinears with light-like separation 
(e.g. [96]). Applying the decomposition in (B.2) to the vector Dirac structure leads to the 
following parameterisation 


{V{p)\qi{0)'y^,q2{z)\0) = flmy Jdu e^^^^ | ^p^, {u) p^gl{u) - g||,3(^) | , 

(B.3) 

for a vector meson V[qiq 2 ] coupling to a light-like separated vector quark-bilinear.^® Consid¬ 
ering all Dirac structures this amounts to a total of eight DAs <?i||,_L,51^3 and 

^®The latter are often denoted by n± or [n, n] with the two remaining directions being labelled by _L. 

2 ®Above the Wilson line between 0 and z, rendering the matrix element gauge invariant, is omitted for brevity. 
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which loosely follow the nomenclature of the nucleon parton distributions functions. It is read¬ 
ily verified, using rjp = 0, f{u) = 1 for / = 4>\\: 9±j 9\\,3 and (B.2), that in the limit z —)■ 0 

the left hand side of (B.3) reduces to rj^fymv as required. Using (B.l) and (B.2) this can be 
written in terms of the actual momentum of the vector meson p and the near light-like distance 
X as follows^^ 


(U(p) 1 91(0)7^92(3;) 1 0) = flmv ~ (^<('||(«)-ffl(^) + ;[^"iy 3 ;^(^||('' 3 )-G)[(u)) 


px my 


+ (gliu) + (5||,3(^^)-29l(n) + ())||(u)) 

1 


2{pxY 

(^(p)l9i(0)cr/,,,92(x)|0) = - ify jdue^^P^ \ {p^py - p^p^) ((/>±(u) -h ^ myX^K^{u)) 


+ {Pf,Xy - pyX^) my ^/i||*^(u) - ^(j)±iu) - 

I 

+ ^ iVf^Xy - PuXfj) (hx^siu) - 0±(u))| , 


(U(p)|9i(0)7^7592(x)|0) ^([(n)-h 


2 2 
mf/X 




16 


Gliu) , 


(U(p)|9i(0)92(x)|0) 


= 7,fvivx)my / , 


(B.4) 


where the notation py 3 = 93 and h ±^3 = /13 has been introduced (wrt to Ref. [96]) in order to 
declare the polarisation of the DAs. Additionally 


= (1 - d+)h\^^ , gliu) = (1 - 6+)9Uu) , G1 = (1 - 5+)Gl , (B.5) 

take into account valence quark mass corrections 

^ ^ 91 ) _ /y (^92 ^ ^gi) 

(/y"ly) (/I'.my) 

consistent with the normalisation 

= 1 , 4> = {4>\\,±,9T^ht’^^^9\\,3,h±,3} , hi^’Ku) = [ dv(f){v), (B.7) 

Jo 

and the EOM of the LCDAs. The twist-4 meson mass corrections Ay _i_ (wrt 4>\\,± DAs) have 
been introduced in [98]. The twist-5 meson mass corrections G][’“ (wrt 9 ^^^ DAs) are intro¬ 
duced in this work for the first time. The DAs G]])’^ and Ay _|_ are not subject to a particular 
normalisation whereas A [Ay —G][](l) = 0 is necessary and partly motivated the reinvestigation 
of the twist-4,5 DA in this work. More details follow just below. We will see in section B.2 
that the four additional DAs Ay _i_ and G][’“ can be obtained from the eight basic ones (at least 
in the asymptotic limit). 

^^Note that (U(p)|gi(0)75q'2(2;)|0) = 0 by parity conservation of QCD. 

^®With due apologies we follow the notation in [96] and not the newer and more systematic notation introduced 
in [97] because of reasons of familiarity. A dictionary between the two notations is shown in Tab. 12. 

The definition of G7 is adapted such that Ay remains as in [98] and (0||, g±) and (Ay, Gx) take on an 
analoguous role. 
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A striking feature are the {px)~^- and (px)“^-ternis which originate from the change of 
variables (B.l) and (B.2). From a conceptual viewpoint the limit px —)• 0 ought to exist and 
imply conditions on the DA which have to be obeyed automatically by the solutions. 

More concretely the 1/px factors are removed by 


1 

px 




=0 


Jo 


(B.8) 


with the normalisation condition as indicated. In Eq. (B.4) this concerns the following five 
functions combintations 


1 

i{px) 

{■/'ll - 9l) 



= h[h - 


I 

i^{px 

)2('T 



= Hhf - 

- - ]^hx,o]{u. 

1 

i{px) 

(/i±,3 - </>±) 



= A [^±,3 

- , 

1 

f‘{px 

^2 (£'11,3 2c?_l + (/)||]) 

- 


II 

- 2g'l_ + , 

1 

i{px) 

(A, - Gl) 


A||^^(n) 

= 7i[A|| - 

-Gl](u), 


where ^ = 2n — 1 and I 2 = Iio Ii is a double application of (B.8). The asymptotic form of the 
DAs has been indicated and the dots stand for non-asymptotic corrections. As stated earlier 
all DAs are such that unwanted boundary terms disappear which is guaranteed provided that 


/i[<).](l) = 0, I2[</.](1) = 0, (B.IO) 

with the first and both conditions applying to the case where the DA in (B.9) is written in 
terms of Ii and I 2 respectively. The integrated DAs are those that appear in the explicit results 
quoted in section A.4.2. 


B.2. Determination of m^-LCDA in asymptotic limit 

Introducing more LCDA means that more information is needed to solve for the DAs. We did 
not systematically aim to do this but present an alternative and possibly new way to determine 
the asymptotic form of the DA directly from (B.3). We expand all quantities systematically 
to first order in rriy, using Eq. (B.l), including in particular the pz in exponential factor in 
(B.3), 

(B.U) 

By matching the first power in my this leads to the following identifications at the level of 
asymptotic DAs (cf. footnote 29) 

A||(u) = - 47i[^((/>||)](u) + 4/2[4ffl-5||,3 - 3(/>||](u) = 24 u^m2 + ... , 

A±(u) = - 47 i[^())_l](u) + 472[/i_l,3 - (/>_l](u) = 24u^m^ + ... , 

G)[(u) = — 47i[^p)[](u) , = 6uu{l — uu) + ... , 
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DA 

4>± 

h 

9l 



hf 

4" 

4^) 

A_i_ 

A|| 

G1 

G1 

NN 

4>2y 



'^iv 


^\y 

V'ly 

^iv 

4>iy 

<i>\y 

<^5,y 

^iv 

form 

- 

- 

3.9 

3.9 

3.7 

3.7 

3.24 

4.21 

4.22 

- 

- 

- 

expl 

(B.14) 

(B.14) 

3.18 

3.17 

3.16 

3.15 

3.27 

4.24 

4.25 

{B.12) 

{B.12) 

{B.12) 

Dirac 

^ py 

7/. 

7 m 

7 m 75 

1 

^ pu 

7 m 

^ pu 

^ py 

7 m 

7 m 

7 m 75 

twist 

2 

2 

3 

3 

3 

3 

4 

4 

4 

4 

5 

5 


Table 12: Translation table between old [96,98] and new notation (NN) [46,97]. The third and fourth 
line indicate the reference to the formal and explicit solution of the DAs. The twist-3,4 
DAs refer to Refs. [46, 97] respectively. The first two DAs are of leading twist 2 and are 
not referenced since they have been known for a long time e.g. [95]. The last three DAs 
are obtained in this work. The second last line denotes the Dirac structure of the DA with 
l,(T^,y and 7/^75 being chiral odd and chiral even respectively. 

Gl(«) = - 4/i[C(7l](n) , = 12^2^2 + ... , (B.12) 

where the dots stand for non-asymptotic corrections. This method, convenient as it is, can 
not determine non-asymptotic corrections since one would need to expand in the z-coordinate 
of the quarks as well as the Wilson line. This leads to higher dimensional local operators and 
3-particle DAs which are both non-asymptotic. Hence if only the asymptotic DAs are required 
then we do not need to do this expansion. We might turn to a more systematic study of 
this method in future work. Our confidence in this alternative method, for determining the 
asymptotic DAs, is borne out of several consistency checks. 

• The DA A|| is such /i[A|| — G][](l) = 0 which guarantees the, previously discussed, finite 
limit p-x —)• 0. The latter allows for the substitution (B.9), i.e. the 1/px-pole is removed 
by partial integration.^^ 

• The asymptotic form A± from the literature is reproduced. 

• We adapt the moment equation [96] [Eq. 4.7] for to G^’’^ which can be done by 
replacing (n -|- 2) —>■ (n -|- 4) where the extra additive factor of 2 originates from the 
extra power of in the LCDA expansion. The adapted moment equation in the limit 
of vanishing quark mass and no 3-particle DA becomes 

^(n + 4)Mf“)=MF), £duCGliu), = 2u - 1) . (B.13) 

It is readily verified that the asymptotic DAs given in (B.12) satisfy the moment equation 
exactly. 


remind the reader that /i [Ay |bbl](1) A 0 was our motivation to investigate the DAs. It should be men¬ 
tioned that it is possible that the EOM can be satisfied to the given order in my by using the ultrarelatvistic 
approximatoin ri{p) —>■ p/m for Ay. For consistency one should use the same approximation to determine 
the light-cone DAs. 
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Last but not least the obtained asymptotic DAs do verify the EOM (5) at the level of 
the correlation function and therefore FFs. This was our original motivation to look into 
this matter. 


B.3. Explicit DAs used for this work 

In this section we provide the actual DAs used in this paper to the given approximation. 
For more complete solutions for the lower DAs we refer to the references in Tab. 12. The 
conventional approach for twist -2 and twist-3 DAs is the expansion in conformal spin (e.g. 
Gegenbauer moments) analogous to the partial wave expansion of SO{3). For the twist 4 there 
is the conformal spin expansion as well as a renormalon model e.g. [97]. In this work we only 
solve for the asymptotic DA for twist 4 which is the lowest order in the conformal expansion. 
We estimate the effect of this to be at the 1%-level which is well beyond the uncertainty. 

Q /o 

For the twist-2 DA we expand up to second order in the Gegenbauer polynomials Cn! {C = 
2u-l) 

4>±,\\{u) = 6u{l-u) + , (B.14) 

which is a standard approximation in view of the lack of reliable knowledge on higher moments. 
The twist-3 DAs are obtained from the twist-2 and the leading twist-3 3-particle DAs [96] 


g±=u dv ^ +u dv —— 
Jo V V 


1 


9± = 


= - 


L4o 


L4o 


V Ju V _ 
dv^+ f'dv^ 


+ 6uu c| + - 4/2)^ 

+ 5 C| ( 3^2 - 1) + ^(cal - ca|/2) (3 - 30^^ + 35^4^ ^ 
+ — 30^^ + 35^^) , 


= u I dv —+ u I dv —+ 6uu ( y 


'^hiv) 


"^hiv) 


-,3/2, 


(B.15) 


with ^g{u) = 2(()||(u) -I- 5+^4''j_{u) + 6-4>'j_{u) and ^hiu) = ‘^(p±iu) - 5+(0'j_(u) - ^/24>'j_{u)) + 
l/26-(f)'j_{u). The contributions of the 2-particle DAs are given implicitly whereas the 3- 

particle DA-contributions have been given explicitly. The four parameters C 3 , , ca| and ca| 

are G-parity even parameters of the three twist-3 3-particle DAs as given in [Eq. 3.11] in [46]. 
SU (3)-breaking parameters can be neglected at the current level of precision. 

As stated above for the twist-4 and twist-5 DAs we employ the asymptotic form which means 
that we set the Gegenbauer moments 01^2 j the 3-particle DA parameters and the quark masses 
to zero with respect to the more general solution. The asymptotic twist-4 DAs are given by 


/i_L^3 = Quu , g'||_3 = 6 un , A_l = 2AvJu^ , Ay = 24n^ri^ (B.16) 


and the newly introduced twist-5 DAs 


G]]_ = 6 uu(l — uu) , = 12 rt^ri^ . 


(B.17) 


The determination of Ay and G][’“ are new and discussed in the previous section. The numerical 
input is given in Tab. 1 except for the values for ci;|“, cul and a)| which are taken from table one 
in [46] and are related to the parameters previously used (e.g. [2]) as follows 3/2(}^(jj^ = ca|, 
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3/2Qa;J = a;|[ and l/2Qa;^ = For the sake of completeness and clarity we give the form 
of the integrated asymptotic DAs given in (B.9): 

g'll^3(u) = — . (B.18) 

C. Decay constants from experiment 

C.l. The Neutral Decay Constants /po ^ from —>• e‘*‘e“ 

We improve the discussion on the extraction of the decay constants of the p^, oj and cj) from 
—>■ e'’'e“ with respect to the earlier work [74], The effects on the decay constants due 
to mixing are taken into account at the level of matrix elements. Previously the mixing was 
abstracted from the state mixing. The relation between the two is commented on in section 

C.1.3. 

The three vector mesons p^,uj and cj) are flavour neutral and can therefore mix into each 
other through QCD and QED. The mixing of (^-o; is driven by QCD, p-ui is due to QED and 
mu,d-quark mass difference whereas (j)-p requires both forces and can therefore be neglected. 
In —)■ e'’“e“ the meson couple to the electromagnetic current as follows 

jT = (c-1) 

with Qu,d,s being the charges of the quarks and the quark currents are defined by 

^ • (C-2) 

The label I stands for isospin as well as ideal mixing (i.e. cj) being a pure ss-state). The 

currents and are of isospin 1 = 0 and 7 = 1 for respectively. 

It is our goal to extract the following decay constants 

, cl{uj\q-if,q\{)) = , cJo(/3°|97/^910) = V^^rUpof'^'’ , (C.3) 

relevant for the description of flavour transition via the weak force. Above p denotes the 
polarisation vectors, the || superscript on fy is omitted and c“o = —c^o = = cf, = \/2 are 

prefactors taking into account the quark composition of the wave functions. The effect of the 
mixing is investigated in a two step procedure of (p-uo and p-u mixing. 

C.1.1. Effective couplings to the electromagnetic current 

mixing: In order to asses the effect of the (j)-u mixing the following matrix elements 
are needed 

. (C.4) 

They have been computed to be e = e{mfj,) ~ e{muj) = 0.05(2) in the pioneering papers of QCD 
sum rules [14]. Note, the effect is driven by contributions of four quark condensates, estimated 
in the vacuum saturation approximation, and we have therefore assigned a conservative 40% 
error to e. 
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p uj mixing: The analogous p, ui decay constants have been computed in reference [99] 
by using finite energy sum rules. Effects are due to different QED and m^^d-quark mass 
differences.^^ Their results, neglecting the p-cj) mixing, translates into the notation analogous 
to (C.4) as follows: = \/6E® ~ 5.9(12) MeV and /£" = V2F^ = -4.8(10) MeV. We have 

enlarged the uncertainty in view of possible duality violation of finite energy sum rules [100]. 


C.1.2. Scaling factors due to mixing 


We parameterise the mixing effects in terms of correction factors, denoted by k, to the matrix 
element to the electromagnetic current, 

, (p°|j7|0) = , 

■ (C-5) 

Using the numerical input of the previous section we get 










( f7Qu + Qd\ 

I 4 ' ) 

(. fi' v^Q. \ 

y Qu + Qd j 

( f^° Qu-Qd \ 

\ f^^Qu + Qd) 

( 1 F Qd I 

( ']^,qFq.) 



Qu T \ 

V2Qs J 

V2Qs \ 
Qu T Qd j 


~ 1/(1.037(16)) , 
~ 1/(0.933(28)) , 
~ 1/(1.08(2)) , 

~ 1/(0.990(2)) . 


(C.6) 


The impact of the mixing on the extraction of the decay constants is heavily affected by the 
charge ratios V^Qsj{Qu + Qd) = —\/2 and (Qu — Qd)/{Qu + Qd) = 3. The reader is reminded 
that p-cj) mixing is neglected since it requires the strong as well as the electromagnetic force 
which is expected to be a small effect. 

The experimental branching ratios are [80] 

BR(p° ^ e+e")=(4.72 ± 0.05) x 10"® , BR(a; ^ e+e") = (7.28 ± 0.14) x 10"^ , 

BR{(j) e+e") =(2.95 ± 0.30) x 10"® . (C.7) 


The theoretical expression for the decay rate is given by 


r(U° ^ e+e") 


47r o? 
3 my 


fyCy + O 



M^J ’ 


(C.8) 


where the coefficients cy in the limit of no mixing (i.e. k —)• 0) can be read-off from (C.l) 
CpO = {Qu - Qd)‘^l‘2' = 1/2, = {Qu + Qd)^l'2 = 1/18 and = Qs = 1/9. The effect of 

^^Note computing the QED corrections to the local matrix element, i.e. which we call decay constant, is not the 
same as computing the QED corrections to the corresponding leptonic decays themselves. The latter are more 
complex, requiring the computation of virtual and real corrections taken into account in the experimental 
analysis. 
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mixing leads to the following replacements 

CpO , Cuij , 

and results in a shift of p , a; and (p decay constant of roughly —1%, 0.5% and 4% respectively. 
It is noticed that the individual effects of the u-p and the uj-cp mixing are around ±8% but 
do almost cancel each other out. 

Including the mixing effects we get the following decay constants for the currents (C.2) 

= (215.6 ± 2Br ± Ir^ ± ± 0 ) MeV = 216(3) MeV, 

= (196.5 ± 2Br ± Ir. ± 4^p ± 6^) MeV = 197(8) MeV, 

4^ = (233.0 ± 2Br ± Ir^ ± 3^ ± 0 ) MeV = 233(4) MeV, (C.9) 

where the uncertainties in the other input parameters are irrelevant. Errors are added in 

quadrature in the final result. To clarify the notation we quote the example ~ 

0.99/pO with /pO from (C.8). 

Finally we get for the p and oj decay constants coupling directly to u and d quark currents 

/p 0 ^=/p 0 + /p7=(215.6 + 5.9) MeV = 221.5(3) MeV , 

4?=/po° - /p7=(215.6 - 5.9) MeV = 209.7(3) MeV , 

=(196.5 - 4.8) MeV = 191.7(8) MeV , 

- /<^J°=(196.5 + 4.8) MeV = 201.3(8) MeV , (C.IO) 

where we have taken the same uncertainties as in (C.9). 

C.1.3. Comment on state mixing versus decay constant mixing 

The mixing of states and decay constants are related but can be quantitatively different. 
The former is one of the effects contributing to the latter. Below we present evidence that in 
reality the mixing of states dominates the mixing of the lo — cp decay constants. 

For example if one assumes that ~ {SU (3)jT’-symmetry for which there is empirical 

evidence), | \ <C | fPtj \ and takes into account the cp-u) state mixing 

\u) ^ \uji) - eui(i>\(pi), \(p) ^ \(pi) + e^^\uji), (C.ll) 

one arrives at 

rrioj , s m0 i n^ 

^ e{m^) , ~ e[muj) , (C-12) 

rriuj 

which is reasonably well satisfied. A recent determination of the mixing angle by the KLOE 
collaboration is given by euxf, = 3.32(9)° ~ 0.58(2). Using, as previously [14], e = e^m^p) ~ 
eirricc) = 0.05(2), Eq. (C.12) is equivalent to 0.45 ~ 0.5(2) and 0.7 ~ 0.5(2) which is satisfied 
within errors. It is to be concluded that the effect of (p-uj decay constant mixing is driven by 
the state-mixing. 

One could put forward the same procedure for the p-uj system but there are complications. 
The p-LO system is more delicate since the closeness of the two states means that the mixing 

^^This is particularly enhanced for rj-rj' system because of the effect of the chiral anomaly. 
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angle is effectively a complex number because diagonal and off-diagonal self energies are com¬ 
plex. The off-diagonal self energy acquires an imaginary part through the isospin violating 
a; —7- vrvr —?■ p transition; a circumstance which has been neglected in the literature for a long 
time! The off-diagonal self energy has been determined to be ~ (—4620 ± 220niodei i 

170data) + (—6100 ± ISOOmodei i 11 lOdata)* MeV^ [101] by using a recent BaBar analysis [102], 
The value of the mixing angle through ep^^ = Yipi^{m?p)/((rriuj — ir^/2)^ — (rup — irp/2)^) then 
comes with a large error; especially on the real part which is decisive. The small error on 
previous determinations turned out to be an artefact of neglecting the imaginary part of the 
off-diagonal self energy [99] . In view of this situation we chose to directly use the computations 
on the mixing of the decay constants and abandon the mixing of state picture. 

C.2. Charged decay constants from r"*" —decays 

The same standard procedure is applied as in [74] including in addition the sizeable leading 
electroweak corrections (due to a ln(m^/mr)-term) [103,104]. Implementing this amounts 
to making the replacement BR(r+ —)• = 1-015 • BR(t''' —)• K*~^u) Using the 

input parameters jUis] — 0.225, jUudj — 0.974, BR(r''' —)• = 25.22 ± 0.33 and BR(r''' —)• 

K*~^v) = 1.20 ± 0.07, a{mv) — 1/135.4 and the r lifetime Tr = (290.3 ± 0.5) • 10“^®s [80] we 
get 

fx* = 204(7) MeV , fp+ = 210(4) MeV . (C.13) 

C.3. Final results summarised 

In view of the many details discussed and numbers quoted we summarise our results for the 
reader’s convenience. The difference in the charged and neutral p decay constants is 6 MeV 
and we therefore choose to average them fp^ = 213(5) MeV slightly enhancing the uncertainty. 

The final results for the decay constants coupling to the currents (C.2) are then taken from 

(C.9), (C.13) and the above mentioned average 

=213(5) MeV, =197(8) MeV, 

/^^=233(4)MeV , /k* =204(7) MeV . (C.14) 

The decay constant for the pure flavour currents to the p^- and w-meson are (C.IO) 

7 ^0^=221.5(3) MeV , =209.7(3) MeV 

/i“)=191.7(8)MeV /i'’')=201.3(8) MeV . (C.15) 

In our tables and computation we will use the decay constants (C.14) omitting the additional 
labels. The results of fx* is consistent with [104], is treated similarly to [105], whereas 
our discussion on p-u mixing is more detailed in terms of explicit results. We would like to 
add a comment concerning QED corrections. The experimental analyses are performed using 
photon showers (e.g. Photos [106]) and subtracting the large part of the final state photons. 
A fully consistent treatment of QED corrections might be carried out in the future for which 
we may expect a global shift (i.e. multiplicative factor in front of all decay constants) at or 
below the 1%-level. 
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As for the flavour specific decay constants we leave it to the reader to scale the B —)• p^,uj 
FFs appropriately. For example for 


1.040 , 

^ 0.973 , (C.16) 

where F stands for any of the seven FFs. The scale factors k are = f^t^/fp and 

~ which upon using (C.14) and (C.15) amount to 

fc(pO,,) = 1.040 , V,<i) = 0-985, /^(^,„) = 0.973 , = 1.022 . (C.17) 

Scaling the FF as in (C.16) is a reasonable procedure since, in practical computations, all 
input parameters of the DA are made dependent on the normalisation of the longitudinal 
decay constant. 

D. Conversion between form factor bases 

D.l. Helicity basis 

In this appendix we give the projection of the FFs onto the helicity basis which is convenient 
for the computation of angular observables. Using the Jacob Wick polarisation tensors^^ we 
define: 


= ep{p){K*{p,ri{m{p)))\sT^^h\B{pB)) , m{t) = m(0) = 0 , m(±) = ± (D.l) 


where p = 0 , ±, f is the polarisation index which is not summed over and F^ = ± 75 ), 

Fy = 7^(1 T 75 ) and F^ = {2i D)^( 1 ± 75 ) correspond to tensor, vector and derivative FFs. We 
get 


X(ll) 

X(0) 

xW 


^(X(+) - a(-)) = iV2y^)Xi , 

4=(A(+) + XF)) = ±iV2{ml - m\*)X2 , 
v 2 


. v^(m| + 3m^ 
2mK* 


2 


Xp, 



(D.2) 


where Xq = X 2 — C 2 'i{q^)X^ with 


C23{q^) 

A(g^) 


_ A(g^) 

(m^ - rn\^.){m\ + 3m|.* 

{{mp + mK*f - 


r) 


- mK*f - q^) , 


(D.3) 


and A being the Kallen-function. We infer that at the kinematic endpoint where A = 0, only 
the X 2 structure contributes in accordance with general findings on endpoint symmetries [77]. 

^®Cf. appendix A [77] where the polarisation tensors 77 and t correspond to 7 and /3 respectively. 
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For X = T, V,T), Xi is given by Tj (with Tp = 0) V* and Vi in Eq. (3) and (10) respectively. 
The relation of Tq and Vo to T 23 and A 12 used in the literature (e.g. [10]) is as follows: 


rp ^_ -16mpm|-, 

° (mp + mp'*)(m| + 3m|.* — g^) ° 


where 


Ai 2 = 


T 23 = ■ 

We further notice that 


(ms + rriK*)'^ {m?B “ "^1'* “ “ Kq‘^)^2 

16msm|.* (ms + rriK*) 

g^/2(m| + 3m|-, - q‘^)Ai + \{q^)mK -/[thb + rnK*)A 3 

8 msm|.*("iB - u^K*) 

(m| - m|.«) (m| + 3m^. - g^) r 2 - X{q‘^)T 3 
SmBm^-t. {niB — niK*) 


Ai 2 ( 0 ) = -^^ 3 ( 0 ) = -^^o( 0 ) , 

oTTlBirax* oTTIb'^K* 

which we implement, besides ri( 0 ) = 22 ( 0 ), into the fit as a constraint. 


(D.4) 


(D.5) 


(D. 6 ) 


cTfiuq’' 



75 

type 

Ti,2,3 

V 

^1,3,0 

^0 

traditional 

Ti,2,3 

Vi 

V2,3,P 

Vp 

EOM (6)-(9) 

^±,11 ~ ^1,2 ,To ~ T 23 

Vy ~ Vi 

V[| ~ V 2 , Vo ~ Ai2 , Vp 

Vp 

helicity 


Table 13: The conversion factors between the traditional and the EOM FFs is given in (3). The 0- 
helicity FF are given by r[V]o( 9 ^) = T[V] 2 ( 9 ^) —C 23 ( 9 ^)T[V] 3 (g^) with the kinematic function 
given as in (D.3) and the dependence of the factor relating r[V]_Lj| ^ T[V]i ,2 can be read 
off from (D.2). The 0 helicity FFs T 23 and A 12 whose notation is inspired by (D.5) are 
related to Tq and Vq as given in (D.4). 


D.2. Overview of form factor notation 

Not including the derivative FFs there are seven independent FFs of which all others are linear 
combinations. The basis Ti^ 2 , 3 , V and ^ 1 , 3,0 is the traditional basis (e.g. [2]; note: A 2 is 
linearly dependent on ^ 0,3 cf. (3)). The basis Ti^ 2 , 3 , and Vi, 2 , 3 ,p is suited for the EOM and 
the conversion between the two is given in (3). The helicity basis T[V]_lj|^o and Vp is suited 
for phenomenology with r[V]j_j| ~ 2"[V’]i,2 and T[V]o = r[V] 2 (Q'^) — C 23 {q'^)T[V] 3 {q'^). The 
0-helicity EFs ^12 and r 23 have been introduced in [11] and their relation to the traditional 
basis is given in (D.5). An overview is given in Tab. 13. 


E. Plots of form factors as a function of 2 : 

The plots of the FEs in the z-variable can be found in Figs. 5,6, and 7 for the modes B —)• K*, 
Bg ^ (j) and Bg —)• K*, respectively. 
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z(q2) z(q2) 

Figure 5: Combined LCSR and lattice fit to i? —>■ K* FFs, where lattice data points are indicated in 
red, LCSR points in blue, the gray solid band shows the combined 3-parameter fit and the 
red dashed band the 2-parameter lattice fit from Ref. [11]. 
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Figure 6: Combined LCSR and lattice fit to Sg —>■ </) FFs, where lattice data points are indicated in 
red, LCSR points in blue, the gray solid band shows the combined 3-parameter fit and the 
red dashed band the 2-parameter lattice fit from Ref. [11]. 
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Figure 7: Combined LCSR and lattice fit to Bg —t K* FFs, where lattice data points are indicated in 
red, LCSR points in blue, the gray solid band shows the combined 3-parameter fit and the 
red dashed band the 2-parameter lattice fit from Ref. [11]. 
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F. SSE coefficients 


In this appendix we list the central values and uncertainties of the SSE expansion coefficients 
of the B —>■ K*, B ^ p, B ^ ui, Bg ^ (j) and Bg —>■ K* FFs from LCSR (Tab. 14) as well 
as the combined fits to LCSR and lattice data for the B —)• K* and Bg ^ cf) FFs (Tab. 15). 
Note that of the 21 parameters for each transition, two are in fact redundant due to the exact 
relations (17).^^ 

In addition to these central values and uncertainties, we also provide the full correlation and 
covariance matrices as ancillary files downloadable from the arXiv preprint page. The data 
are contained in 5 JSON files named [Process] _ [Fit] . json, where [Process] is BKstar for 
B —^ K*, Brho for B ^ p, Bomega for R —)• cj, Bsphi for Bg ^ cp and BsKstar for Bg —)• K* 
FFs; [Fit] is LCSR for the fit to LCSR only (valid at low q'^) and LCSR-Lattice for the 
combined fit valid in the full range. 

The JSON format can be easily used in Mathematica. For example, reading in the file for 
the B-^ K* LCSR FFs, 

data = Import ["BKstar_LCSR.json"] 

the central value of can be accessed simply via 

OptionValue[data, "central" -> "Tl" -> "aO"] 

and the correlation between and as 

OptionValue[data, "correlation" -> "AOV" -> "ala2"] 

and similarly for the objects "uncertainty" and "covariance". In Python, the corresponding 
commands would read 

import json 

with open(’BKstar_LCSR.jsonO as file: 
data = json.load(file) 

and the parameters can be accessed via 

data['central’][’Tl’][’aO’] 
data[’correlation’][’AOV’][’ala2’] 

etc. 


G. Lifetime effect in Bg — 

To compare the experimental measurement of the Bg —)• cppp branching ratio and angular 
observables from an untagged data sample to the theoretical predictions, the difference in 
finite width AT^ between the Bg mass eigenstates of widths Tl and Th has to be taken into 
account [84,107]. This leads to a difference between experimentally accessible time-integrated 


^"^Due to this redundancy, the 21 x 21 covariance matrices do not have full rank. Invertible covariance matrices 
can be obtained by removing the two redundant rows and columns. 
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B - 


K* 

B 

—) 

■ P 

B 

—> 

UJ 

Bs 

— 

^ </> 

Bs 


K* 

“o 

0.36 

± 

0.05 

0.36 

± 

0.04 

0.33 

± 

0.05 

0.39 

± 

0.05 

0.31 

± 

0.05 


-1.04 4 

I 0.27 

-0.83 4 

:0.20 

-0.83 4 

I 0.30 

-0.78 4 

I 0.26 

-0.66 4 

:0.23 

Aq 

1.12 

± 

1.35 

1.33 

± 

1.05 

1.42 

± 

1.25 

2.41 

± 

1.48 

2.57 

± 

1.44 

Ai 

«0 

0.27 

± 

0.03 

0.26 

± 

0.03 

0.24 

± 

0.03 

0.30 

± 

0.03 

0.23 

± 

0.03 


0.30 

± 

0.19 

0.39 

± 

0.14 

0.34 

± 

0.24 

0.48 

± 

0.19 

0.27 

± 

0.19 

Ai 

"2 

-0.11 4 

I 0.48 

0.16 

± 

0.41 

0.09 

± 

0.57 

0.29 

± 

0.65 

0.13 

± 

0.56 

Ai2 

«o 

0.26 

± 

0.03 

0.30 

± 

0.03 

0.27 

± 

0.04 

0.25 

± 

0.03 

0.23 

± 

0.03 


0.60 

± 

0.20 

0.76 

± 

0.20 

0.66 

± 

0.26 

0.76 

± 

0.20 

0.60 

± 

0.21 


0.12 

± 

0.84 

0.46 

± 

0.76 

0.28 

± 

0.98 

0.71 

± 

0.96 

0.54 

± 

1.12 


0.34 

± 

0.04 

0.33 

± 

0.03 

0.30 

± 

0.04 

0.39 

± 

0.03 

0.30 

± 

0.03 


-1.05 4 

I 0.24 

-0.86 4 

:0.18 

-0.83 4 

I 0.29 

-1.03 4 

I 0.25 

-0.90 4 

:0.27 

«2^ 

2.37 

± 

1.39 

1.80 

± 

0.97 

1.72 

± 

1.24 

3.50 

± 

1.55 

2.65 

± 

1.33 


0.28 

± 

0.03 

0.27 

± 

0.03 

0.25 

± 

0.03 

0.31 

± 

0.03 

0.24 

± 

0.02 


-0.81 

14 

I 0.19 

-0.74 4 

:0.14 

-0.72 4 

I 0.22 

-0.87 4 

I 0.19 

-0.76 4 

:0.20 

af 

1.95 

± 

1.10 

1.45 

± 

0.77 

1.41 

± 

1.01 

2.75 

± 

1.19 

2.08 

± 

1.00 


0.28 

± 

0.03 

0.27 

± 

0.03 

0.25 

± 

0.03 

0.31 

± 

0.03 

0.24 

± 

0.02 

ap 

0.40 

± 

0.18 

0.47 

± 

0.13 

0.41 

± 

0.23 

0.58 

± 

0.19 

0.34 

± 

0.19 

a? 

0.36 

± 

0.51 

0.58 

± 

0.46 

0.46 

± 

0.57 

0.89 

± 

0.71 

0.52 

± 

0.61 

723 

«o 

0.67 

± 

0.08 

0.75 

± 

0.08 

0.68 

± 

0.09 

0.68 

± 

0.07 

0.60 

± 

0.08 

To 3 

1.48 

± 

0.49 

1.90 

± 

0.43 

1.65 

± 

0.62 

2.11 

± 

0.46 

1.58 

± 

0.56 

To 3 

Ot2 

1.92 

± 

1.96 

2.93 

± 

1.81 

2.47 

± 

2.19 

4.94 

± 

2.25 

3.65 

± 

3.27 


Table 14: Fit results for the SSE expansion coefficients in the fit to the LCSR computation only. These 
numbers are provided (to higher accuracy) in electronic form along with the full correlation 
matrices as arXiv ancillary files. The LCSR FFs are usually taken to be valid in the range 
from 0 to 14GeV^. 


CP-averaged observables Oexp and the theoretical definition of CP-averaged observables Otheo 
in the flavour eigenstate basis. The former and the latter are defined as 

Oexp= / —0(t), Otheo = 0(t = 0) , (G.l) 

Jo TBs 

where tBs is the lifetime of the Bg 

0{t)(Bs —t IT ) ~ 2 4’IT^IT ) + (J{Bs{t) —7- (j)ii^li )] . (G.2) 

In the case where only vector operators are present (i.e. i^eff &7At(75)'S^7'"(75)^), the time- 
dependent CP-averaged observables 0{t) can be written as functions Fq of bilinears of time- 
dependent transversity amplitudes Jbx,aY{t) 

0{t) = Fo{JbX,aY{t)) , (G.3) 


55 




B - 


K* 

Bs 

— 

f <P 

Bs 


K* 

< 4 ° 

0.37 

± 

0.03 

0.42 

± 

0.02 

0.36 

± 

0.02 

of” 

-1.37 d 

I 0.26 

-0.98 d 

I 0.24 

-0.36 d 

:0.20 

Ao 

“2 

0.13 

± 

1.63 

3.27 

± 

1.36 

8.03 

± 

2.07 


0.30 

± 

0.03 

0.29 

± 

0.01 

0.22 

± 

0.01 


0.39 

± 

0.19 

0.35 

± 

0.10 

0.24 

± 

0.16 

Ai 

“2 

1.19 

± 

1.03 

1.70 

± 

0.79 

1.77 

± 

0.85 

Ai2 

“o 

0.27 

± 

0.02 

0.27 

± 

0.02 

0.27 

± 

0.02 


0.53 

± 

0.13 

0.95 

± 

0.13 

1.12 

± 

0.11 


0.48 

± 

0.66 

2.15 

± 

0.48 

3.43 

± 

0.78 

< 

0.38 

± 

0.03 

0.36 

± 

0.01 

0.28 

± 

0.02 

ar 

-1.17 d 

I 0.26 

-1.22 d 

I 0.16 

-0.82 d 

:0.19 

ar 

2.42 

± 

1.53 

3.74 

± 

1.73 

5.08 

± 

1.42 


0.31 

± 

0.03 

0.30 

± 

0.01 

0.24 

± 

0.01 


-1.01 d 

I 0.19 

-1.10 d 

I 0.08 

-0.75 d 

:0.15 


1.53 

± 

1.64 

0.58 

± 

1.00 

2.49 

± 

1.37 


0.31 

± 

0.03 

0.30 

± 

0.01 

0.24 

± 

0.01 

a? 

0.50 

± 

0.17 

0.40 

± 

0.08 

0.31 

± 

0.15 

a? 

1.61 

± 

0.80 

1.04 

± 

0.61 

1.58 

± 

0.93 

To 3 

«o 

0.67 

± 

0.06 

0.65 

± 

0.04 

0.60 

± 

0.04 

To 3 

1.32 

± 

0.22 

2.10 

± 

0.33 

2.40 

± 

0.27 

To 3 
^2 

3.82 

± 

2.20 

6.74 

± 

1.80 

9.64 

± 

2.03 


Table 15: Fit results for the SSE expansion coefficients in the combined LCSR + lattice fit. These 
numbers are provided (to higher accuracy) in electronic form along with the full correlation 
matrices as arXiv ancillary hies. 


Jtx.aY(t) = Ai{t)Al{tr + Ai{t)Al{tr . (G.4) 

where a, 6 = 0, ||,T are the vector meson polarisation indices and X,Y = L,R denote the 
chirality structure of the lepton production. The CP-conjugated amplitude is 

-</>^) , (G.5) 


where 7?|| q = +1 and r]± = —1 are the CP-eigenvalues of the amplitudes and 
refers to the conjugation of all weak (CP-odd) phases. Defining the quantity 


cA _ _ —i(ps " "Q 

“’a ^ ’ 

where (ps is the Bg mixing phase, one can write 
JbX,aY{t) = JbX,aY{0)-^ 




(G.6) 


(G.7) 
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with 


,bX,aY _ 

•^Ar — 




In summary for the experimental and theoretical expression in (G.l), we obtain 


(G.8) 


JbX,aY\exp 

JbX,aY\theo 


dt 


Jo TB, 
JbX,aYiO) , 


JbX,aY{t) — 


1 + Us-J^ 


bX,aY 

Ar 




-JbX,aY{0) , 


(G.9) 


where 


_ hi - Th 
^^~Tl + Th 


Ar, 

2r, 


1 

TBs 


-{Tl + Th) . 


(G.IO) 


As a simple example, we consider the differential branching ratio at low in the SM, within 
naive factorisation and the heavy quark and massless lepton limit, where the transversity 
amplitudes read 


= \/2NmB^{l - s) 



2mbmB ^eff^ 

A 





(G.ll) 


with s = /m?B ^ being a normalisation factor including the CKM elements VtbV^g. Note 
that the soft FFs _i_ are not to be confused with the ratio of amplitudes in (G.6). One finds 


= = ^^’-^^ = +1. (G.12) 

The theoretical and experimental GP-averaged differential branching ratios in the assumed 
limit read 


dBR 


dq^ 


= TB, 


theo 


+ |A|f p + |A|f p + 


= rB. [2{\Ai\^ + \Af\^) + \A^\^ + \A§\^] 


(G.13) 


and 


dBR 




a=±,||,0 ^ 


= BBs 


1 




I _ y2 M---L1 ■ 1---LI / ■ 


(G.14) 


where 

A a aclL.clL AdR.aR 

Aaf = -^Ar = “^Ar = -da ■ 

At low the sizeable longitudinal polarization fraction of the (/>-meson signals that the last 
term in (G.14) dominates, so the time-integrated branching ratio at low q^ is suppressed by 
0{ys), where Us = 0.62(5) [81], with respect to the prompt one. This is in agreement with the 
findings in [84]. 
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